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Dyson-Schwinger equations furnish a Poincare covariant framework within which to 
study hadrons. A particular feature is the existence of a nonperturbative, symmetry 
preserving truncation that enables the proof of exact results. The gap equation reveals 
that dynamical chiral symmetry breaking is tied to the long-range behaviour of the 
strong interaction, which is thereby constrained by observables, and the pion is precisely 
understood, and seen to exist simultaneously as a Goldstone mode and a bound state of 
strongly dressed quarks. The systematic error associated with the simplest truncation 
has been quantified, and it underpins a one-parameter model efficacious in describing an 
extensive body of mesonic phenomena. Incipient applications to baryons have brought 
successes and encountered challenges familiar from early studies of mesons, and promise 
a covariant field theory upon which to base an understanding of contemporary large 
momentum transfer data. 
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1. Introduction 

A central goal of contemporary nuclear physics is to understand the properties 
of hadrons in terms of the elementary excitations in quantum chromodynamics 
(QCD): quarks, gluons and ghosts, because this is a crucial step in validating the 
theory. Here "elementary excitations" means those quantum fields in terms of which 
QCD's Lagrangian is naturally expressed; for example, the analogues in quantum 
electrodynamics (QED) are the electron and photon. Just as positronium does not 
appear in the QED Lagrangian but is understood as a relativistic bound state in 
the theory, hadrons do not appear in QCD's Lagrangian: they are supposed to be 
bound states of quarks and gluons. It is an observational fact that hitherto only 
colourless hadrons have been observed and that hadrons cannot be "ionised;" i.e., 
unlike positronium, for which the addition of a few (~ 7) electron-volts will separate 
the electron and positron, no amount of energy available presently is sufficient to 
break a hadron into separated coloured constituents. To understand this unique 
feature of confinement is one of the most significant challenges in physics. 

The defining problems of hadron physics can now be stated: specify and solve 
the quantum field theoretical bound state problem whose solution is the hadron 
spectrum; calculate the interactions of these hadrons between themselves and with 
electroweak probes; and use these interactions as incisive tools with which to identify 
the origins of confinement and elucidate its effects on, and expression in, observables. 
These problems are essentially nonperturbative. 

In considering hadron physics it is natural to think of numerical simulations of 
lattice-QCD, which over the last thirty years has become an independent branch 
of high-energy physics with numerous large-scale collaborations. An essential as- 
pect of lattice-QCD is the enumerable, finite grid used to represent the spacetime 
continuum, which means that a comparison between simulations and experiment 
necessarily involves extrapolations. In modern studies the extrapolation to the con- 
tinuum limit; i.e., to small spacing between the lattice sites, is reliably handled. 
However, that is not true of the extrapolation to infinite volume. Furthermore, 
most simulations continue to employ the so-called quenched approximation, which 
introduces an a priori unquantifiable systematic error. Significant resources are be- 
ing expended to overcome these difficulties, by developing improved algorithms, and 
by the obvious expedient of acquiring faster computers with more memory. Lattice- 
QCD's current status can be appraised from the proceedings of any of the major 
lattice conferences, e.g., Refs. pQ. There are naturally successes and problems, and 
we leave their description to practitioners, e.g., Ref. [2]. 

In describing many aspects of hadron physics, light-front concepts and tech- 
niques are useful, as may be seen, e.g., from Ref. The application, exploration 
and improvement of these methods, too, forms an identifiable subfield but one that 
is not insular: much is being gained by capitalising on the opportunities that ex- 
ist to inform other nonperturbative approaches and models, and the feedback that 
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brings. This framework has seen numerous recent successes, e.g., a renormalisation 
group approach to the light-front HamiltonianElapplicationJHl of light-front Hamil- 
tonian techniques to the transverse lattice formulation of QCD, and refinements 
of discretised light-front quantisation.^ A perspective on the approach's status and 
future is presented in Ref. 

Reviewing this material it is readily apparent that drawing a connection be- 
tween QCD and hadron observables is difficult, and that is why modelling remains 
a keystone of hadron physics. Constituent-quark- like models currently provide a 
peerless description of the baryon spectrum and baryon decaysP^ correlating a 
wealth of data via few parameters. Furthermore, while such models do not simul- 
taneously give a satisfactory description of low- lying light-quark mesons, for which 
an accurate representation of dynamical chiral symmetry breaking (DCSB) is es- 
sential, they can be ad apted and applied to guide the search for exotic and hybrid 
hadrons" and glueballs ! ' ' * When focusing on properties of the baryon octet 
and decuplet, such as elastic and transition form factors at small to moderate mo- 
mentum transfer, mean field models of baryon structure ameliorate some of the 
deficiencies of constituent-quark models and provide additional insig 
At higher momentum transfers a proper expression of Poincare covariance becomes 
important in models of such processes! I I 

Contemporary Dyson-Schwinger equation (DSE) studies complement these ap- 
proaches. Modern comparisons with and predictions of experimental data can prop- 
erly be said to rest on model assumptions but they can be tested within the frame- 
work and also via comparison with lattice-QCD simulations, and the predictions 
are very good. Furthermore, progress in understanding the intimate connection 
between symmetries and truncation schemes has enabled the proof of exact results. 
Herein we review recent phenomenological applications and the foundation of their 
success. There are naturally challenges, which we shall also highlight. 



2. Dyson-Schwinger Equations 

The best known DSE is the simplest: the Dyson or gap equation, which describes 
how the propagation of a fermion is modified by its interactions with the medium 
being traversed. In QCD that equation is: h 

S- 1 (p)^Z 2 (i 1 -p + m harc )+ Z 1 J^ g 2 D^(p-q)^S{q)T a „(q;p); (2.1) 

i.e., the renormalised DSE for the dressed-quark propagator. In Eq. I|2.1() . D^(k) 
is the renormalised dressed-gluon propagator, T®(q;p) is the renormalised dressed- 
quark-gluon vertex, rribarc is the A-dependent current-quark bare mass that appears 
in the Lagrangian and J := J d q/iJZir) 4 represents a translationally-invariant 

a ln fact, these quantum mechanical models must be adapted because exotic hadrons are, by defi- 
nition, those states which cannot be constructed in the simplest versions of the model, with only 
constituent-quark degrees of freedom. 

'We employ a Euclidean metric throughout, with: {7 M , 7^} = 2<5 M „; 7^ =7^; and a-b = ^^t—i a i^i' 
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regularisation of the integral, with A the regularisation mass-scale. c In addi- 
tion, Zi(( 2 , A 2 ) and ^(C 2 ; A 2 ) are the quark-gluon- vertex and quark wave function 
renormalisation constants, which depend on the renormalisation point, £, and the 
regularisation mass-scale, as does the mass renormalisation constant 

Z m (C 2 , A 2 ) = Z 4 (C 2 ,A 2 )/Z 2 (C 2 , A 2 ). (2.2) 

The renormaliscd mass is m(£) :— rrib are (A) / Z m (( 2 , A 2 ). When ( is very large the 
right-hand-side (r.h.s.) can be evaluated in perturbation theory whereby one finds 

vn 

with 7 m = 12/(33 — 2iV/), where Nf is the number of current-quark flavours that 
contribute actively to the running coupling, and Aqcd and rh are renormalisation 
group invariants: rh is the renormalisation-group-invariant current-quark mass. 
The solution of Eq. 1)2.1(1 is the dressed-quark propagator, which takes the form 

S-\p) = z 7 • P A(p 2 , C 2 ) + B(p 2 , C 2 ) = I [i 7 • p + M{p 2 )] (2.4) 

and is obtained by solving the gap equation subject to the renormalisation condition 
that at some large spacelike ( 2 

S-HP)\ P 2 =C , =*7-P + m(C). (2.5) 

The gap equation illustrates the features and flaws of each DSE. It is a non- 
linear integral equation for S(p) and hence can yield much-needed nonperturbative 
information. However, the kernel involves the two-point function D^Qc) and the 
three-point function T®(q;p). The equation is therefore coupled to the DSEs these 
functions satisfy. Those equations in turn involve higher n-point functions and 
hence the DSEs are a tower of coupled integral equations with a tractable problem 
obtained only once a truncation scheme is specified. It is unsurprising that the best 
known truncation scheme is the weak coupling expansion, which reproduces every 
diagram in perturbation theory. This scheme is systematic and valuable in the anal- 
ysis of large momentum transfer phenomena because QCD is asymptotically free 
but it precludes any possibility of obtaining nonperturbative information, which we 
identified as a key feature of the DSEs. 

In spite of the problem with truncation, gap equations have long been used effec- 
tively in obtaining nonperturbative information about many-body systems as, e.g., 
in the Nambu-Gorkov formalism for superconductivity!^ The positive outcomes 
have been achieved through a simple expedient of employing the most rudimentary 
truncation, e.g., Hartree or Hartree-Fock, and comparing the results with obser- 
vations. Naturally, agreement under these circumstances is not an unambiguous 

c It is only with a translationally invariant regularisation scheme that Ward-Takahashi identities can 
be preserved, something that is crucial to ensuring vector and axial-vector current conservation. 
The final stage of any calculation is to take the limit A — » oo. 
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indication that the contributions omitted are small nor that the model expressed in 
the truncation is sound. However, it does justify further study, and an accumulation 
of good results is grounds for a concerted attempt to substantiate a reinterpretation 
of the truncation as the first term in a systematic and reliable approximation. 

2.1. Nonperturbative truncation 

To explain why Eq. (|2.1|l is called a gap equation we consider the chiral limit, which 
is readily defined^ because QCD exhibits asymptotic freedom and implemented by 
employing^ 1 

Z 2 (C 2 ,A 2 )m baro (A) = 0, A»C- (2.6) 

An equivalent statement is that one obtains the chiral limit when the renormali- 
sation-point-invariant current-quark mass vanishes; i.e., rh = 0. In this case the 
theory is chirally symmetric, and a perturbative evaluation of the dressed-quark 
propagator from Eq. (|2.1|) gives 



b° pc Ap 2 ) 





\p 2 ' 

9 






m 





■°0; (2.7) 



viz., the perturbative mass function is identically zero in the chiral limit. It follows 
that there is no gap between the top level in the quark's filled negative-energy Dirac 
sea and the lowest positive energy level. 

However, suppose that one had at hand a truncation scheme other than per- 
turbation theory and that subject to this scheme Eq. I|2.1|l possessed a chiral limit 
solution B°(p 2 ) ^ 0. Then interactions between the quark and the virtual quanta 
populating the ground state would have nonperturbatively generated a mass gap. 
The appearance of such a gap breaks the theory's chiral symmetry. This shows that 
the gap equation can be an important tool for studying DCSB, and it has long been 
used to explore this phenomenon in both QED and 

The gap equation's kernel is formed from a product of the dressed-gluon prop- 
agator and dressed-quark-gluon vertex. However, in proposing and developing a 
truncation scheme it is insufficient to focus only on this kernel. We have observed 
that the gap equation can be a tool for studying DCSB but it is only useful if the 
truncation itself does not destroy the chiral symmetry. 

Chiral symmetry is expressed via the axial-vector Ward-Takahashi identity: 

P f ,r 5ti {k;P) = S-\k + )i l5 + i 75 S- 1 (k_), k± = k±P/2, (2.8) 

wherein T^^^k; P) is the dressed axial- vector vertex. This three-point function sat- 
isfies an inhomogeneous Bethe-Salpeter equation (BSE): 

P)] tu = Z 2 [ 75 7mL + [ A [S(q + )T^(q; P)S(q-)UKZ(q, k; P) , (2.9) 



in which K(q,k] P) is the fully- amputated quark-antiquark scattering kernel, and 
the colour-, Dirac- and flavour-matrix structure of the elements in the equation is 
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denoted by the indices r,s,t,u. The Ward-Takahashi identity, Eq. 1)2. 8[l . entails 
that an intimate relation exists between the kernel in the gap equation and that in 
the BSE. (This is another example of the coupling between DSEs.) Therefore an 
understanding of chiral symmetry and its dynamical breaking can only be obtained 
with a truncation scheme that preserves this relation, and hence guarantees Eq. 1)2.8(1 
without a fine-tuning of model-dependent parameters. 



2.1.1. Rainbow-ladder truncation 

At least one such scheme exists. 1 ^ Its leading-order term is the so-called renor- 
malisation-group-improved rainbow-ladder truncation, whose analogue in the many 
body problem is an Hartree-Fock truncation of the one-body (Dyson) equation com- 
bined with a consistent ladder-truncation of the related two-body (Bethe-Salpeter) 
equation. To understand the origin of this leading-order term observe that the 
dressed-ladder truncation of the quark-antiquark scattering kernel is expressed in 
Eq. ijZHJl via 



r BM («;P)] tt 'f := [S(q + )T 5ll (q;P)S(q-)} sr K™(q,k;P) 



= -g 2 (C 2 )D P Ak-q) [T a p (k + ,q + )S(q + )] u , [S(q.)T a a (q.,k.)] u , u [T 5fi (q; P)] Vu , 

(2.10) 

wherein we have only made explicit the renormalisation point dependence of the 
coupling. One can exploit multiplicative renormalisability and asymptotic freedom 
to establish that on the kinematic domain for which Q 2 := (k — q) 2 ~ k 2 
large and spacelike 



q 2 is 



[L(q,k;P)]^' =-4ira(Q 2 )D^ c (Q) 



^7p5 frcc (g + ) 



(2.11) 

where a(Q 2 ) is the strong running coupling and, e.g., S 66 is the free quark prop- 
agator. It follows that on this kinematic domain the r.h.s. of Eq. 1)2.11(1 describes the 
leading contribution to the complete quark-antiquark scattering kernel, K^(q, k; P), 
with all other contributions suppressed by at least one additional power of l/Q 2 . 
The renormalisation-group-improved ladder-truncation supposes that 



K™{q,k-P) - -Aira(Q 2 )D s ^{Q) 









Y 7p_ 


ts 





(2.12) 



is also a good approximation on the infrared domain and is thus an assumption 
about the long-range (Q 2 < 1 GeV 2 ) behaviour of the interaction. Combining Eq. 
((2.12(1 with the requirement that Eq. 1(2.8(1 be automatically satisfied leads to the 
renormalisation-group-improved rainbow-truncation of the gap equation: 

S-\p) = Z 2 (i 1 -p + m barc ) + ^a{Q 2 )Df v c {p - q)^- 7 „ S(q) ^ . (2.13) 

As will become apparent, the rainbow-ladder truncation provides the foundation 
for an explanation of a wide range of hadronic phenomena. There are some notable 
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exceptions; e.g., the description of the scalar meson sector is unconvincing, and it 
is natural to seek the reason. The recognition tha t th is truncation is the first term 
in a systematic procedure has provided an answer P2I 

2.2. Systematic procedure 

The truncation scheme of Ref . [221 is a dressed-loop expansion of the dressed-quark- 
gluon vertices that appear in the half-amputated dressed-quark-antiquark scattering 
matrix: S 2 K, a renormalisation-group invariant!^ All n-point functions involved 
thereafter in connecting two particular quark-gluon vertices are fully dressed. The 
effect of this truncation in the gap equation, Eq. (|2.1[1 . is realised through the 
following representation of the dressed-quark-gluon vertex, iY a ^ = ^X a T p = ^T^: 

+ -f I 9 2 D^(£) D T , T (l + k-p) 7 r< S(p - t) la , r^(£, -k, k-p) + [...}. 

(2.14) 

Here T 3g is the dressed-three-gluon vertex and it is readily apparent that the lowest 
order contribution to each term written explicitly is 0(<? 2 ). The ellipsis repre- 
sents terms whose leading contribution is 0(g 4 ): viz., the crossed-box and two-rung 
dressed-gluon ladder diagrams, and also terms of higher leading-order. 

This expansion of S 2 K, with its implications for other n-point functions, yields 
an ordered truncation of the DSEs that guarantees, term- by-term, the preserva- 
tion of v ector and axial-vector Ward-Takahashi identities, a feature that has been 
ex ploitecPE3EHl 

to prove Goldstone's theorem and other exact results in QCD. It 
is readily seen that inserting Eq. 1)2.14(1 into Eq. (|2.1|l provides the rule by which 
the rainbow-ladder truncation can be systematically improved. 

2.2.1. Planar vertex 

The effect of the complete vertex in Eq. (I2.14f) on the solutions of the gap equation 
is unknown. However, insights have been drawn from a study^of a more modest 
problem obtained by retaining only the sum of dressed-gluon ladders; i.e., the vertex 
depicted in Fig. 12.11 The elucidation is particularly transparent when one employs 
the following choice for the dressed-gluon line in the figured" 

2V(*) := 9 2 D^(k) = (Spv - (2ir) 4 Q 2 S A (k) , (2.15) 

which defines an ultraviolet finite model so that the regularisation mass-scale can 
be removed to infinity and the renormalisation constants set equal to one. d This 



d The constant Q sets the model's mass-scale and using Q = 1 simply means that all mass- 
dimensioned quantities are measured in units of Q. 
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Fig. 2.1. Integral equation for a planar dressed-quark-gluon vertex obtained by neglecting contri- 
butions associated with explicit gluon self-interactions. Solid circles indicate fully dressed propa- 
gators. The vertices are not dressed. (Adapted from Ref. 1231 .) 

model has many positive features in common with the class of renormalisation- 
group-improved rainbow-ladder models and its particular momentum-dependence 
works to advantage in reducing integral equations to algebraic equations with similar 
qualitative features. There is naturally a drawback: the simple momentum depen- 
dence also leads to some model-dependent artefacts, but they are easily identified 
and hence not cause for concern. 

The general form of the dressed-quark gluon vertex involves twelve distinct scalar 
form factors but using Eq. (|2.15l) only IV (p) := IV(p,p) contributes to the gap 
equation. This considerably simplifies the analysis since, in general, 

IVG 5 ) = ai(p 2 )7 M + a 2 {p 2 )l -PPh - a 3 (p 2 ) ipp, + a^p 2 ) 17^7 -p. (2.16) 

The summation depicted in Fig. 12. II is expressed via 

r P (p) = 7^ + 1 i P s( P ) r„(p) s( P ) lp , (2.17) 

and inserting Eq. Q2.1fi|l into Eq. Ij2.17|l one finds ot± = 0, and hence the solution 
simplifies: 

IV (p) = ai(p 2 )7 M + a 2 {p 2 )j-pp ll - a 3 (p 2 )ip M . (2.18) 

The three surviving fun ctions are those which are most important in the dressed- 
quark-photon vertexP3 

One can re-express this vertex as 

00 00 

r M (p) = T » = E H(p 2 ) > + "2b 2 ) 7 • pp, - 4(p 2 )i Pli ] , (2.19) 

1=0 8=0 

where the superscript enumerates the order of the iterate: I 1 * is the bare vertex, 

a\ = 1 , a° 2 = = a° 3 ; (2.20) 

r^ =1 is the result of inserting this into the r.h.s. of Eq l|2.17|l to obtain the one-rung 
dressed-gluon correction; T l = 2 is the result of inserting r^ =1 , and is therefore the 
two-rung dressed-gluon correction; etc. A key observation^ is that each iterate is 
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Fig. 2.2. Recursion relation for the iterates in the fully-resummed dressed-gluon-ladder vertex: 
filled circles denote a fully-dressed propagator or vertex. Using Eq. 12.151 . p+ = p~. (Adapted 
from Rcf. 1551 .1 



related to its precursor via the simple recursion relation depicted in Fig. 12.21 and, 
substituting Eq. I|2.19|l . this recursion yields (s = p 2 ) 




a ^( 8 ):=\ al +1 ( S ) = 0(s: A, B) a % (s) , (2.21) 



I 1 I A ° ° 

0(s;A,B) = --^\ 2A 2 sA 2 — B 2 2AB |. (2.22) 



4AB 4s AB 2(B 2 - sA 2 ) 
A = sA 2 (s) + B 2 (s). It follows that 



a 



= {E oi )<* = T±o a ° (2 ' 23) 

and hence, using Eq. I|2.2()(l . 



4A 

ai = 



1 + 4A ' 

-8 A 2 1 + 2 A 

012 ~ l + 2(B 2 -sA 2 )-8A 2 1 + 4A' ( ' ' 

-8AB 

a 3 = 



1 + 2{B 2 - sA 2 )-8A 2 



The recursion relation thus leads to a closed form for the gluon-ladder-dressed 
quark-gluon vertex in Fig. 12. II viz., Eqs. I|2.18|l . I|2.24|l . Its momentum-dependence 
is determined by that of the dressed-quark propagator, which is obtained by solving 
the gap equation, itself constructed with this vertex. Using Eq. I|2.15|l . that gap 
equation is 

S-\ P ) = t 7 • p + m + 7 „ S{p) r» (2.25) 
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-4 -2 2 4 

S 

Fig. 2.3. A(s), M(s) obtained from Eqs. l2~24l . i2~26l . 12371 with m = 0.023: solid line. All 
dimensioned quantities are expressed in units of Q in Eq. 12.151 . For comparison, the results 
obtained with the zeroth-order vertex: dotted line; and the one-loop vertex: dashed line, are also 
plotted. (Adapted from Ref. 1251 .) 

and substituting Eq. (|2.18|) gives 

A(s) = 1 + 1 B2 [A (2at - sa 2 ) -Ba 3 ], (2.26) 

B(s) = m+ * n9 [B (4ai + sa 2 ) ~ sA a 3 ] ■ (2.27) 
sA z + B £ 

Equations Ij2.26|l . Ij2.27|l . completed using Eqs. I|2.24|l . form a closed algebraic sys- 
tem. It can easily be solved numerically, and that yields simultaneously the complete 
gluon-ladder-dressed vertex and the propagator for a quark fully dressed via gluons 
coupling through this nonperturbative vertex. Furthermore, it is apparent that in 
the chiral limit, m — 0, a realisation of chiral symmetry in the Wigner-Weyl mode, 
which is expressed via the B = solution of the gap equation, is always admissible. 
This is the solution anticipated in Eq. I|2.7|l . 

The chiral limit gap equation also admits a Nambu-Goldstone mode solution 
whose p 2 ~ properties are unambiguously related to those of the m ^ solution, 
a feature also evident in QCDpSI A complete solution of Eq. (|2.25(l is available 
numerically, and results for the dressed-quark propagator and gluon-ladder-dressed 
vertex are depicted in Figs. 12.31 12. 41 It is readily seen that the complete resummation 
of dressed-gluon ladders gives a dressed-quark propagator that is little different from 
that obtained with the one-loop-corrected vertex; and there is no material difference 
from the result obtained using the zeroth-order vertex. Similar observations apply 
to the vertex itself. Of course, there is a qualitative difference between the zeroth- 
order vertex and the one-loop-corrected result: 012.3 7^ in the latter case. However, 
once that effect is seeded, the higher-loop corrections do little. The scale of these 
modest effects can be quantified by a comparison between the values of M(s = 0) = 













Bare 




1 i 00 p 


\ M(s) 


Recursive 
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0.4 



0.2 







a,(S) 






1 loop 




2 loop 




Recursive 


<z 3 (S) 




«2<S) 





-4 -2 2 4 

S 

Fig. 2.4. cti, i = 1, 2, 3, calculated from Eqs. <2~24l . <2~26l . i2~27l with m = 0.023. These functions 
calculated at one-loop (dotted line) and two-loop (dashed line) are also plotted for comparison. 
(Adapted from Ref. 



B(0)/A(0) calculated using vertices dressed at different orders: 
Ei=o,jvr; ^ = N = l N = 2 N=oc 



M(0) 



1 



1.105 1.115 1.117 



(2.28) 



The rainbow truncation of the gap equation is accurate to within 12% and adding 
just one gluon ladder gives 1% accuracy. It is now important to couple this with an 
understanding of how the vertex resummation affects the Bethe-Salpeter kernel. 



2.2.2. Vertex- consistent Bethe-Salpeter kernel 

The renormalised homogeneous BSE for the quark-antiquark channel denoted by 
M can be expressed 

[T M (k;P)}tu = f [ X M{q;P)]sr[K{k,q;P)]Z (2.29) 

Jq 

where: Tm(&; P) is the meson's Bethe-Salpeter amplitude, k is the relative momen- 
tum of the quark-antiquark pair, P is their total momentum; and 

X M(k)P) = S{k + )T M {k;P)S{k-). (2.30) 

Equation (|2.29|) . depicted in Fig. 12. 51 describes the residue at a pole in the solution of 
an inhomogeneous BSE; e.g., the lowest mass pole solution of Eq. (|2.9|l is identified 
with the pion. (NB. The normalisation of a Bethe-Salpeter amplitude is fixed by 
requiring that the bound state contribute with unit residue to the fully-amputated 
quark-antiquark scattering amplitude: M = K + K(SS)K+ . . ., see, e.g., Ref. |28|.) 

On p. [SI we observed that the automatic preservation of Ward-Takahashi identi- 
ties in those channels related to strong interaction observables requires a conspiracy 
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Fig. 2.5. Homogeneous BSE, Eq. 12.291 . Filled circles: dressed propagators or vertices; K is the 
dressed-quark-antiquark scattering kernel. A systematic truncation of S 2 K is the key to preserving 
Ward-Takahashi identities!! 2 ^ (Adapted from Ref. I23I .1 



between the dressed-quark-gluon vertex and the Bethe-Salpeter kernelP2E31 A sys- 
tematic procedure for building that kernel followJ^I from the observatiorP^I that 
the gap equation can be expressed via 



6T[S] 



= 0. 



(2.31) 



where T[S] is a Cornwall- Jackiw-Tomboulis-like effective action. The Bethe-Salpeter 
kernel is then obtained via an additional functional derivative: 



K' 



SS rs 



(2.32) 



With the recursive vertex depicted in Fig. 12.11 the n-th order contribution to 
the kernel is obtained from the n-loop contribution to the self energy: 



= - / A 2V(P - q) l°7„ S(q)l a r:(q,p). 

Jq 



(2.33) 



Since T fJi (p,q) itself depends on S then Eq. Ij2.32|l yields the Bethe-Salpeter kernel 
as a sum of two terms and hence Eq. (|2.29|) assumes the form 



T M {k;P)= f V^(k-<l)l a J» 

Jq 

where we have used the mnemonic 

oo 

A^(g,fc;P)=^AX(9,fc;P). 



XM (q; P) l a T v {q-,k-) + S(q+) K a Mu (q, fc; P) , 

(2.34) 



(2.35) 



n=0 



While T>^ v also depends on S because of quark vacuum polarisation diagrams, 
the additional term arising from the derivative of T>^ v does not contribute to the 
Bethe-Salpeter kernel for flavour nonsinglet systems, which are our current focus, 
and hence is neglected for simplicity. 

Equation l|2.34|l is depicted in Fig. 12.61 The first term is instantly available once 
one has an explicit form for T™ and the second term, identified by the shaded box 
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Fig. 2.6. Upper panel: BSE, Eq. 12.341 . which is valid whenever T„ can be obtained via a recursion 
relation. Lower panel: Recursion relation for , Eq. 12.361 . (Adapted from Ref. 1231 .') 



in Fig. 12.61 can be obtained via an inhomogeneous recursion relatio: 



+ / V pa (k- q )l b lp S(q + +£-k)l a r:- 1 (q + +£-k, q+ ) XM (q;P)l b la 



This equation is also depicted in Fig. 12.61 Combining the two figures, it is apparent 
that to form the Bethe-Salpeter kernel the free gluon line is attached to the upper 
dressed-quark line. It follows that the first term in Eq. I|2.36|l invariably generates 
crossed gluon lines; i.e., nonplanar contributions to the kernel. The character of the 
vertex-consistent Bethe-Salpeter kernel is now clear: it consists of countably many 
contributions, a subclass of which are crossed-ladder diagrams and hence nonplanar. 
Only the rainbow gap equation, obtained with i = in Eq. Ij2.19|l . yields a planar 
vertex-consistent Bethe-Salpeter kernel, namely the ladder kernel of Eq. (|2.12|) . In 
this case alone is the number of diagrams in the dressed- vertex and kernel identical. 
Otherwise there are always more terms in the kernel. 

2.2.3. Solutions for the it- and p-mesons 

Section 2.2.2 recapitulates on a general procedure that provides the vertex-consistent 
channel-projected Bethe-Salpeter kernel once T™ and the propagator functions: A, 
B, are known. That kernel must be constructed independently for each channel be- 
cause, e.g., h- a Mv depends on Xm(<Z; P)- As with the study of the vertex in Sec. 2.2.1, 
an elucidation of the resulting BSEs' features is simplified by using the model of Eq. 





,q-+k-£) S(q- +k — £) Z V 



+ 



(2.36) 
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Table 2.1. Calculated it and p meson masses, in GeV, quoted with Q = 0.48 GeV, in which case 
m = 0.023 5 = 11 MeV. n is the number of dressed-gluon rungs retained in the planar vertex, see 
Fig. l2.l| and hence the order of the vertex-consistent Bethe-Salpeter kernel: the rapid convergence 
of the kernel is apparent from the tabulated results. (Adapted from Ref. 1231 .') 









M£ =0 


M£ =1 


M%= 2 






m 


= 














T, 


m 


= 0.011 


0.152 


0.152 


0.152 


0.152 


P, 


m 


= 


0.678 


0.745 


0.754 


0.754 


P, 


m 


= 0.011 


0.695 


0.762 


0.770 


0.770 



1(2.15(1 . for then the Bethe-Salpeter kernels are finite matrices [cf. (1 — O) -1 in Eq. 
(|2.23(l ] and the homogeneous BSEs are merely linear, coupled algebraic equations. 

Reference [231 describes the solution of the coupled gap and Bethe-Salpeter 
equations for the tt- and p-mesons in detail. Herein we focus on the results, which 
are summarised in Table 12.11 It is evident that, irrespective of the order of the 
truncation; viz., the number of dressed gluon rungs in the quark-gluon vertex, 
the pion is massless in the chiral limit. This is in spite of the fact that the pion 
is composed of heavy dressed-quarks, as is evident in the calculated scale of the 
dynamically generated dressed-quark mass function: see Fig. 12.31 M(0) « Q w 
0.5 GeV. These observations emphasise that the masslessness of the 7r is a model- 
independent consequence of consistency between the Bethe-Salpeter kernel and the 
kernel in the gap equation. Furthermore, the bulk of the p-ir mass splitting is 
present for m = and with the simplest (n = 0; i.e., rainbow-ladder) kernel, 
which demonstrates that this mass difference is driven by the DCSB mechanism: 
it is not the result of a carefully tuned chromo-hyperfine interaction. Finally, the 
quantitative effect of improving on the rainbow-ladder truncation; i.e., including 
more dressed-gluon rungs in the gap equation's kernel and consistently improving 
the kernel in the Bethe-Salpeter equation, is a 10% correction to the vector meson 
mass. Simply including the first correction (viz., retaining the first two diagrams in 
Fig. 12.1(1 yields a vector meson mass that differs from the fully resummed result by 
< 1%. The rainbow-ladder truncation is clearly accurate in these channels. 

2.2.4. Summary 

It is now evident that a Ward-Takahashi identity preserving Bethe-Salpeter ker- 
nel can always be calculated explicitly from a dressed-quark-gluon vertex whose 
diagrammatic content is enumerable. 6 Furthermore, in all but the simplest case, 
namely, the rainbow-ladder truncation, that kernel is nonplanar. 

While we described results obtained with a rudimentary interaction model in 
order to make the construction transparent, the procedure is completely general. 
However, the algebraic simplicity of the analysis is naturally peculiar to the model. 
With a more realistic interaction, the gap and vertex equations yield a system of 



e That is not true if one employs an Ansatz for the dressed-quark-gluon vertex whose diagrammatic 
content cannot be made explicit; e.g., that of Ref. 1301 . 
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twelve coupled integral equations. The Bethe-Salpeter kernel for any given channel 
then follows as the solution of a determined integral equation. 

We have reviewed those points in the construction of Refs. |22I23| that bear upon 
the fidelity of the rainbow-ladder truncation in the gap equation, and in the vector 
and flavour non-singlet pseudoscalar channels. The error is small. In modelling it 
is therefore justified to fit one's parameters to physical observables at this level in 
these channels and then make predictions for other phenomena involving vector and 
pseudoscalar bound states in the expectation they will be reliable. That approach 
has been successful, as we shall illustrate. 

By identifying the rainbow-ladder truncation as the lowest order in a systematic 
scheme, the procedure also provides a means of anticipating the channels in which 
that truncation must fail. The scalar mesons are an example. Paramctrisations 
of the rainbow-ladder tr uncation, fitt ed as described above, give masses for scalar 
mesons that arc too high!^11 3^3J^4: | ^hat wag thought to be a problem. However, 
we now know this had to happen because cancellations that occur between higher or- 
der terms in the pseudoscalar and vector channels, thereby reducing the magnitude 
of corrections, do not occur in the scalar channel!^ wherein the full kernel contains 
additional attraction!^ Indeed, it follows that an interaction model employed in 
rainbow-ladder truncation that simultaneously provides a good description of scalar 
and pseudoscalar mesons must contain spurious degrees of freedom. Quantitative 
studies of the effect of the higher-order terms have begunP^ESEHEIIl 

The placement of the rainbow-ladder truncation as the first term in a procedure 
that can systematically be improved explains clearly why this truncation has been 
successful, the boundaries of its success, why it has failed outside these boundaries, 
and why straightening out the failures will not undermine the successes. 

2.3. Selected model-independent results 

In the hadron spectrum the pion is identified as both a Goldstone mode, associated 
with DCSB, and a bound state composed of constituent u- and d-quarks, whose 
effective mass is ~ 350 MeV. Naturally, in quantum mechanics, one can fabricate a 
potential that yields a bound state whose mass is much less than the sum of the 
constituents' masses. However, that requires fine tuning and, without additional 
fine tuning, such models predict properties for spin- and/or isospin-flip relatives 
of the pion which conflict with experiment. A correct resolution of this apparent 
dichotomy is one of the fundamental challenges to establishing QCD as the theory 
underlying strong interaction physics, and the DSEs provide an ideal framework 
within which to achieve that end, as we now explain following the proof of Rcf. |19j . 

2.3.1. A proof of Goldstone's theorem 

Consider the BSE of Eq. (|2.29|) expressed for the isovector pseudoscalar channel: 




(2.37) 
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with xiili P) — S{q+)T : ! K {q] P)S(q-) obvious from Eq. H2.30j) and j labelling isospin, 
of which the solution has the general form 



Ti(k;P) 



= t j 75 



iE n (k ] P)+ 1 -PF n (k;P) 



1 -kk-PG^{k\P) + a lxv k^P v H v {k; P) 



(2.38) 



It is again apparent that the dressed-quark propagator, the solution of Eq. (|2.1|) . is 
an important part of the BSE's kernel. 

In studying the pion it is crucial to understand chiral symmetry, and its explicit 
and dynamical breaking. These features arc expressed in the axial- vector Ward- 
Takahashi identity, Eq. (|2.8|l . which involves the axial- vector vertex: 

r-A 



n,(k;P) 



= Z 2 



757^" 



{xUq;P)]srK™( q ,k;P), 



(2.39) 



that has the general form 



n,(k;P) 



75 [iMk; P) + 7 ■ kk^Gnik; P) - k v H R (k; P)} 



2 

+rUk;P) + 



P, 



P 2 + ml 



<{k-P) 



(2.40) 



where F R , G R , H R and P 5/i are regular as P 2 -> -to 2 , P^i^k; P) - 0(P 2 ) and 
4^{k\ P) has the structure depicted in (|2.38(1 . This form admits the possibility of at 
least one pole term in the axial- vector vertex but does not require it. 

Substituting (|2.40|1 into (|2.39|1 and equating putative pole terms, it is clear that, 
if present, (j? (k; P) satisfies Eq. H2.37(l . Since this is an eigenvalue problem that only 
admits a / solution for P 2 = —to 2 , it follows that <^(A;;P) is nonzero only 



for P 2 = —to 2 and the pole mass is mj 



Hence, if K supports such a bound 



state, the axial- vector vertex contains a pion-pole contribution whose residue, r A , 
is not fixed by these arguments; i.e., Eq. I|2.4(J|I becomes 



Ti(fc;P) = -l 5 h f ,F R (k;P)+- r kk^G R (k;P)-a^KH R (k;P)} 



n,(k;P) 



r A Pn 



P 2 



-Ti(k-P). 



(2.41) 



Consider now the chiral limit axial- vector Ward-Takahashi identity, Eq. (|2.8() . 
If one assumes to 2 = in Eq. (|2.41() . substitutes it into the l.h.s. of Eq. 12.8|) along 
with Eq. (|2.4|) on the right, and equates terms of order (P„)° and P„, one obtains 
the chiral-limit relations^ 



r A E„(k;Q) 
F R {k;0) + 2r A F n (k;0) 
G R {k;Q) + 2r A G n [k;0) 
H R {t 1 0) + 2r A H n (t 1 0) 



P(fc 2 ), 
A(fc 2 ), 
2A'{k 2 ) 
0. 



(2.42) 
(2.43) 
(2.44) 
(2.45) 
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We know B(k 2 ) = in the chiral limit [Eq. l(2~7jl ] and that a B(k 2 ) ^ solution of 
Eq. 1|2.1[) in the chiral limit signals DCSB. Indeed, in this case 



M(p 2 



large-p 2 2lT 2 'y rl 



(-m°) 



p 2 In 



p 2 /A 2 



QCD 



(2.46) 



where (qq)° is the renormalisation-point-independent vacuum quark condensate!^ 
Furthermore, there is at least one nonperturbative DSE truncation scheme that 
preserves the axial-vector Ward-Takahashi identity, order by order. Hence Eqs. 
(|2.42() - (|2.45(l are exact quark- level Goldberger-Treiman relations, which state that 
when chiral symmetry is dynamically broken^l 

1. the homogeneous isovector pseudoscalar BSE has a massless, P 2 = 0, 
solution; 

2. the Bethc-Salpcter amplitude for the massless bound state has a term 
proportional to 75 alone, with E^(k; 0) completely determined by the 
scalar part of the quark self energy, in addition to other pseudoscalar 
Dirac structures, F n , G w and H 7Tl that are nonzero; 

3. and the axial- vector vertex is dominated by the pion pole for P 2 ~ 0. 
The converse is also true. Hence DCSB is a sufficient and necessary condition for 
the appearance of a massless pseudoscalar bound state (of what can be very-massive 
constituents) that dominates the axial-vector vertex for P 2 ~ 0. 

2.3.2. A mass formula 

When chiral symmetry is explicitly broken the axial- vector Ward-Takahashi identity 
becomes: 

P^Uk;P) = 5- 1 (fc + )* 75 y +*75y5- 1 (fc_) - 2i m(£) rj(fc; P), (2.47) 



where the pseudoscalar vertex is given by 



H(k;P) 



= Z 4 



75y 



+ 



X i(q;P) K™( q> k;P). (2.48) 



As argued in connection with Eq. 12.39f) . the solution of Eq. I|2.48|l has the form 



iTi(k;P) 



-75 [iE£(k; P)+yPF£ + > V kk-PGZ(k;P) 



+ 0-^ k„P v H%(k;P)] 



rp 



P 2 



2 . 



-r£(*;P), 



(2.49) 



where E^, F£, G^ and are regular as P 2 — > — m 2 ; i.e., the isovector pseu- 
doscalar vertex also receives a contribution from the pion pole. In this case equating 
pole terms in the Ward-Takahashi identity, Eq. lt2~47jl . entails!^ 



r A m 2 = 2m(C) r P {Q . 



(2.50) 
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This is another exact relation in QCD. Now it is important to determine the residues 
Ta an d rp. 

A consideration of the renormalised axial- vector vacuum polarisation shows 

r A S ij P» = U5 ij P» = Z 2 tr [ ir i 757lJl S(q + )Ti(q;P)S(q-); (2.51) 

Jq 

i.e., the residue of the pion pole in the axial- vector vertex is the pion decay constant. 
The factor of Zi on the r.h.s. in Eq. (|2.51|l is crucial: it ensures the result is gauge 
invariant, and cutoff and renormalisation-point independent. Equation (|2.51|l is 
the exact expression in quantum field theory for the pseudovector projection of the 
pion's wave function on the origin in configuration space. 

A close inspection of Eq. (|2.48(l . following its re-expression in terms of the 
renormalised, fully-amputated quark- antiquark scattering amplitude: M — K + 
K(SS)K + ..., yieldsPI 

i^r P = Z 4 tr / A lT l l5 S(q + )TUq;P)S(q-) , (2.52) 



wherein the dependence of Z4 on the gauge parameter, the regularisation mass-scale 
and the renormalisation point is exactly that required to ensure: 1) rp is finite in 
the limit A — > 00; 2) rp is gauge-parameter independent; and 3) the renormali- 
sation point dependence of rp is just such as to guarantee the r.h.s. of Eq. (|2.5U|) 
is renormalisation point independent. Equation (|2.52|) expresses the pseudoscalar 
projection of the pion's wave function on the origin in configuration space. 

Let us focus for a moment on the chiral limit behaviour of Eq. I|2.52ll whereat, 
using Eqs. (232) -(233), one readily finds 

,A 

- (m)c = UrUC) = Zi(C, A) N c tr D / S, h=0 (q) . (2.53) 

Jq 

Equation (|2.53|) is unique as the expression for the chiral limit vacuum quark con- 
densate. It is C-dependent but independent of the gauge parameter and the regu- 
larisation mass-scale, and Eq. I|2.53|l thus proves that the chiral-limit residue of the 
pion pole in the pseudoscalar vertex is (— (qq)®)/ fn- Now Eqs. (|2.50l) . i|2.53[) yield 

(/°) 2 ml = -2m(C) (qq)° C + 0(m 2 ) , (2.54) 

where is the chiral limit value from Eq. (|2.51|) . Hence what is commonly known 
as the Gell-Mann-Oakes-Renner relation is a corollary of Eq. (|2.5U|) . 

One can now understand the results in Table 12.11 a massless bound state of 
massive constituents is a necessary consequence of DCSB and will emerge in any 
few-body approach to QCD that employs a systematic truncation scheme which 
preserves the Ward-Takahashi identities. 

We stress that Eqs. (I2.50() - (|2.52|I are valid for a ny values of the current-quark 
masses and the generalisation to Nf quark flavours i fp0 | ^4 | 2 5| 

f 2 H m 2 H = -(qq) c M C H , (2.55) 
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A4? = +m^ 2 is the sum of the current-quark masses of the meson's constituents; 

,A 

fnP tl = Z 2 tr ±(T H ) T Tolf >S(q + )T H (q;P)S(q-), (2.56) 

J a 



with S — diag(5' u , Sd, S s , . . .), T H a flavour matrix specifying the meson's quark 
content, e.g., = |(A X + i\ 2 ), {A 4 } are iV/-fiavour generalisations of the Gell- 
Mann matrices, and 

,A 

(qq)? = if H Z 4 tr ±(T H ) T l5 S(q + )T H (q; P) S(q-) . (2.57) 

Jq 

Owing to its chiral limit behaviour, (qq)^ has been called an in-hadron condensate. 
In the heavy-quark limit, Eq. 1)2.56(1 yields the model-independent resulJ^EHl 



(2.58) 



i.e., it reproduces a well-known consequence of heavy-quark symmetry!^ A similar 
analysis of Eq. (|2.57|l gives a new result 

- (qq)? = constant + O ( — ) for — . (2.59) 

15 \m H J m H 



Combining Eqs. i|235)l . P35|) . one find J 24 l 25 l 



mH oc rhf for — ~0, (2.60) 

mf 

where rhf is the renormalisation-group-invariant current-quark mass of the flavour- 
nonsinglet pseudoscalar meson's heaviest constituent. This is the result one would 
have anticipated from constituent-quark models but here we have reviewed a direct 
proof in QCD. Equation 12.50fl is a single formula that unifies aspects of light- and 
heavy-quark physics and, as we shall illustrate, can be used to gain an insightful 
understanding of modern lattice simulations. 

3. Foundation for a Description of Mesons 

The renormalisation-group-improved rainbow-ladder truncation has long been em- 
ployed to study light mesons and [Sec. 2] it can be a quantitatively reliable tool 
for vector and flavour nonsinglet pseudoscalar mesons. In connection with Eqs. 
1(2.12(1 . 1(2.13(1 we argued that the truncation preserves the ultraviolet behaviour of 
the quark-antiquark scattering kernel in QCD but requires an assumption about 
that kernel in the infrared; viz., on the domain Q 2 < lGeV 2 , which corresponds 
to length-scales > 0.2 fm. The calculation of this behaviour is a primary challenge 
in contemporary hadron physics and there How- 
ever, at present the efficacious approach is to model the kernel in the infrared, which 
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enables quantitative comparisons with experiments that can be used to inform the- 
oretical analyses. 

The most extensively applied model is specified by usinJ^OE] 



g(Q 2 ) _ 4tt 2 Q 2 /u) 2 | 8^ 2 7 „ 

Q 2 ~ u F ' 



6 



hi 



1 + QVA 2 QCD ) 



T W), (3.1) 



in Eqs. (l2~T2l . (I2~T31) . Here, T[Q 2 ) = [1 - exp(-Q 2 /[4m 2 ])]/Q 2 , m t = 0.5 GeV; 



t = c 2 - 1; 7 m = 12/25; anJ23 A QCD = = 0.234 GeV/ The true parame- 

ters in Eq. I|3.1|) are D and cj, which together determine the integrated infrared 
strength of the rainbow-ladder kernel; i.e., the so-called interaction tension^] a A _ 
However, we emphasise that they are not independent. °^ in fitting to a selection 
of observables, a change in one is compensated by altering the other; e.g., on the 
domain u> £ [ 0.3. 0.5] GeV, the fitted observables are approximately constant along 
the trajectory-**! 

ujD = (0.72 GeV) 3 . (3.2) 

This correlation: a reduction in D compensating an increase in lu, acts to keep a 
fixed value of the interaction tension. Equation l|3.1|l is thus a one-parameter model. 

3.1. Rainbow gap equation 

Inserting Eq. (f 3 . 1 1) into Eq. 1)2.13(1 provides a model for QCD's gap equation and 
in applications to hadron physics one is naturally interested in the nonperturbative 
DCSB solution. A familiar property of gap equations is that they only support 
such a solution if the interaction tension exceeds some critical value. In the present 
case that value ■Ma* -2.5 GeV/fm. This amount of infrared strength is suffi- 
cient to generate a nonzero va cuu m quark condensate but only just. An acceptable 
description of hadrons require J^a A ~ 25 GeV/fm and that is obtained witrl^J 

D = (0.96 GeV) 2 . (3.3) 

This value of the model's infrared mass-scale parameter and the two current- 
quark masses 

m„(lGeV) = 5.5MeV, m s (lGeV) = 125MeV, (3.4) 
defined using the one-loop expression 

to evolve m u (19GeV) = 3.7MeV, m s (19GeV) = 85MeV, were obtained in Ref. 
|51j by requiring a least-squares fit to the n- and K-meson observables listed in 

■^NB. Eq. I5TI gives ct(m? z ) = 0.126. Comparison with a modern valued 0.117 ± 0.002, means 
that a smaller Aqcd is acceptable in the model, if one wants to avoid overestimating the coupling 
in the ultraviolet, but not a larger value. 
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Table 3.1. Comparison of experimental values with results for n and K observables calculated using 
the renormalisation-group-improved rainbow-ladder interaction specified by Eq. 13.11 . quoted in 
MeV. The model's sole parameter and the current-quark masses were varied to obtain these results. 
The best fit parameter values are given in Eqs. 13.31 . 13.41 . (Adapted from Ref. 1511 . 1 





m v 


m K 
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Sk 


CalcP 
Expt.EH 


138 
138 


497 
496 


93 
92 


109 
113 



Table 13.11 The procedure was straightforward: the rainbow gap equation [Eqs. 
Q2.5|) . H2.13fl . (|3.1|) ] was solved with a given parameter set and the output used to 
complete the kernels in the homogeneous ladder BSEs for the tt- and if -mesons [Eqs. 
(j2~T2|) . (|2~T3|) . l(2~37|) . jlSHt with r j for the tt channel and T j T K+ = ±(A 4 + iX 5 ) 
for the Kj. These BSEs were solved to obtain the tt- and if -meson masses, and the 
Bethe-Salpeter amplitudes. Combining this information gives the leptonic decay 
constants via Eq. (|2.56|l . This was repeated as necessary to arrive at the results 
in Table 13.11 which were judged satisfactory. The model gives a vacuum quark 
condensate 

-<M>?GeV = (0-242 GeV) 3 , (3.6) 

calculated from Eq. H2.53JI and evolved using the one- loop expression in Eq. (|3.5|l . 

With the model's single parameter fixed, and the dressed-quark propagator ob- 
tained, it is straightforward to compose and solve the homogeneous BSE for vector 
mesons and calculate properties analogous to those in Table l3~T1 The prediction:]^ 
of the model are reproduced in Table 13.21 The expression in QCD for a vector 
meson's electroweak decay constant ia^l 

fl Ml = \Z 2 tr [ A (T H r 1 ,S( q+ )r*( q ;P)S( q -), (3.7) 

Jq 

where Ml is the meson's mass. This quantity characterises decays such as p — ► 
e + e~, t — » K*v Tl and in the heavy-quark humP"! 

flac-Lp, fl*f H , (3.8) 

V M H 

reproducing additional familiar consequences of heavy-quark symmetry. 

Given the discussion in Sec. 2.2, the phenomenological success of the rainbow- 
ladder kernel, evident in the results in Tables 13.11 and 13.21 is unsurprising and, 
indeed, was to be expected. 

3.2. Comparison with lattice simulations 

Since the solution of the gap equation has long been of interest in grappling with 
DCSB in QCD, in Figs. 13. ll 13.21 we depict the scalar functions characterising the 
renormalised dressed-quark propagator: the wave function renormalisation, Z(p 2 ), 
and mass function, M(p 2 ), obtained by solving Eq. I|2.13|l using Eq. I|3.1I) . The in- 
frared suppression of Z(p 2 ) and enhancement of M(p 2 ) are longstanding predictions 
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Table 3.2. Experiment cf. predictions of a rainbow-ladder kernel for simple vector meson observ- 
ables. No parameters were varied to obtain the results. The root-mean-square error over predicted 
quantities is < 8%. NB. A charged particle normalisation is used for /X in Eq. 13,71 . which differs 
from that in Eq. 12.561 by a multiplicative factor of y/2. (Adapted from Ref. 1511 . ) 
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742 
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1072 
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Expt. 53 


771 


892 


1019 


217 


227 
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Rel. -Error 


0.04 


-0.05 


-0.05 


0.05 
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of DSE studies'^ which could have been anticipated from Ref. |55|. This prediction 
has recently been confirmed in numerical simulations of quenched lattice-QCD, as 
is evident in the figures. 

It is not yet possible to reliably determine the behaviour of lattice Schwinger 
functions for current-quark masses that are a realistic approximation to those of 
the u- and d-quarks. Therefore a lattice estimate of m^, f^, (qq)° is absent. To 
obtain such an estimate, Ref. |49j used the rainbow kernel described herein and 
varied {D,uj) in order to reproduce the lattice data. A best fit was obtained with 

D = (0.74 GeV) 2 , to = 0.3 GeV , (3.9) 

at a current-quark mass of 0.6m] GoV Ri 14ra u [Eq. (|3.4|l ] chosen to coincide with 
that employed in the lattice simulation. Constructing and solving the homogeneous 
BSE for a pion-like bound state composed of quarks with this current-mass yields 

m m,~14m u = q 48 GeV ^ ym a ~14m„ = Q.094 GeV . (3.10) 

The parameters in Eq. (|3.9|l give chiral limit results@^ 

f° = 0.068 GeV, -(qq)° 1GcW = (0.19 GeV) 3 , (3.11) 

whereas Eqs. Q3.2[l . I|3.3|l give /° = 0.088 GeV, which agrees with the estimate of 
chiral perturbation theory, and the vacuum quark condensate in Eq. I|3.6|) . Following 
this preliminary study it is important to reanalyse the lattice data using different 
models for the infrared behaviour of the scattering kernel, K 7 so that the quenching 
error can be estimated. 

The unification of light and heavy pseudoscalar meson masses via the mass for- 
mula in Eq. (|2.50|l has also been quantitatively explored using the rainbow-ladder 
kernel. That is illustrated in Fig. 13.31 wherein the calculated mass of a uq pseu- 
doscalar meson is plotted as a function of m q (£), with m u (() fixed via Eq. (|3.4I) . 
The calculations are depicted in the figure by the solid curve, which is, in MeVj^-l 

m H = 83 + 500VX + 310 X, X = m c q / A QCD . (3.12) 

The curvature appears slight in the figure but that is misleading: the nonlinear term 
in Eq. (|3.12|) accounts for almost all of m w (the Gell-Mann-Oakes-Renner relation 
is nearly exact for the pion) and 80% of rriK- NB. The dashed line in Fig. 13.31 
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Fig. 3.1. Wave function renormalisation. Solid line: solution of the gap equation using Eqs. 
13. II . 13.31 : data: lattice simulations)*^! obtained with m = 0.036/a ~ 60MeV; dashed-line: gap 
equation solution using Eqs. 13.11 . 13.91 . The DSE study used a renormalisation point £ = 19 GeV 
and a current-quark mass 0.6 ml GcV [Eq. <3~4T so as to enable a direct comparison with the lattice 
data. (Adapted from Ref. 1491 .1 



fits the K, D, B subset of the data exactly. It is drawn to illustrate how easily 
one can be misled. Without careful calculation one might infer from this apparent 
agreement that the large-m g limit of Eq. 1)2.50)1 is already manifest at the s-quark 
mass whereas, in reality, the linear term only becomes dominant for m q > 1 GeV, 
providing 50% of mry and 67% of tub- The model predicts mJ GeV = 1.1 GeV, 
m\ GcW = 4.2 GeV, values typical of Poincare covariant treatmentsPSI 

Equation 1)3.12)1 can be used as a basis for extrapolating the results of lattice 
simulations to realistic values of the light current-quark masses. It is a true ex- 
pression of essential consequences of DCSB and the importance of incorporating 
such constraints in the analysis of lattice data is beginning to be appreciated!^ 
The rainbow-ladder kernel has also been employed in an analysis of a trajectory 
of fictitious pseudoscalar mesons, all composed of equally massive constituents!^ 
(The only physical state on this trajectory is the pion. ) The DSE study predicts!^ 

Mm - 2m " =2.2, (3.13) 

m H m=ms 

in agreement with a result of recent quenched lattice simulations!^ It provides an 
understanding of this result, showing that the persistent dominance by the term 
nonlinear in the current-quark mass o wes itself to a large value of the in-meson 
condensates for light-quark mesons; e.g.p^ (qq)i g c v = (—0.32 GeV) 3 , and thereby 
confirms the large-magnitude condensate version of chiral perturbation theory. (NB. 
This last observation is also supported by Eq. I|3.11|l and the associated discussion.) 
References |57I58| provide vector meson trajectories too. 
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Fig. 3.2. Mass function. Solid line: solution of the gap equation using Eas. 13.11 . 13.31 : data: lattice 
simulations|^l obtained with m = 0.036/a ~ 60MeV; dashed-line: gap equation solution using 
Eqs. 13.11 . 13.91 . The DSE study used a renormalisation point ( = 19 GeV and a current-quark 
mass 0.6mJ GcV [Eq. J-j"""}]. (Adapted from Ref. 14*51 .1 



3.3. Pion's valence-quark distribution function 



The momentum-dependent dressing of quark and gluon propagators is a fact. It is 
certainly the keystone of DCSB, materially influences hadron observables and quite 
likely plays a central role in confinement. This was anticipated in the Global Colour 
Model !^ l 93 l v 4 l As we shall subsequently illustrate, the rainbow-ladder kernel is 
unique today in providing a direct description and unification of a wide range of me- 
son phenomena in terms of a single parameter that characterises the long-range be- 
haviour of the quark-antiquark interaction. Establishing this is a reward for substan- 
tial effort. However, it is a recent development. Historically, algebraic parametri- 
sations of the dressed-quark propagator and Bethe-Salpeter amplitudes, based on 
the known behaviour of numerical solutions, were used to e xpedite a comparison be- 
tween theory and experiment. Successes with that approach 
provided the impetus for refining the DSE method and its direct numerical appli- 
cations. 

The utility of an algebraic form for the dressed-quark propagator and Bethe- 
Salpeter amplitudes is self-evident: calculating the simplest elastic form factor re- 
quires the repeated evaluation of a multidimensional integral whose integrand is a 
complex- valued function, and a functional of the propagator and amplitudes. The 
expedient remains of use, notably in connection with baryons, and can be illus- 
trated b y r eviewing a recent calculation of the pion's valence-quark distribution 
function."^ 

The cross section for deep inelastic lepton-hadron scattering can be interpreted 
in terms of the momentum-fraction probability distributions of quarks and gluons 
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Fig. 3.3. Solid line: pseudoscalar uq meson's mass as a function of rriq, f = 19GeV, with a fixed 
value of m£ corresponding to m^ GeV = 5.5 MeV, Eq. 13.41 . The experimental data points are 
from Rcf. 1531 as are the errors assigned to the associated heavy-quark masses. The dashed curve 
is a straight line drawn through the K, D, B masses. (Adapted from Ref. 1561 . See also Ref. 1581 .) 



in the hadronic target, and since the pion is a two-body bound state with only u- 
and d-valence-quarks it is the least complicated system for which these distribution 
functions can be calculated. However, in the absence of pion targets, their measure- 
ment is not straightforward and they have primarily been inferred from Drell-Yan 
measurements in pion-nucleus collisions!^ 

The distribution functions provide a measure of the pion's quark-gluon sub- 
structure but they cannot be calculated perturbatively. Fortunately, DSEs furnish 
a sound theoretical description of pion properties, providing a model- independent 
explanation of its essentially dichotomous nature as both a Goldstone boson and a 
low-mass bound state of massive constituents, and also supply an efficacious phe- 
nomenological tool. Hence they are ideal for exploring the nature of the pion's 
parton distribution functions. The valence-quark distribution function, u^(x), is 
of particular interest because its pointwise behaviour is affected by aspects of con- 
finement dynamics; e.g., by the mechanisms responsible for the finite extent and 
essentially nonpointlikc nature of the pion. (NB. u 71 = in the C/-parity sym- 
metric limit of QCD.) 

3.3.1. Handbag contributions 

Reference [7S] focuses on a calculation of the "handbag diagrams" illustrated in Fig. 
I5~l which are the only impulse approximation contributions to virtual photon-pion 
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Fig. 3.4. "Handbag" contributions to the virtual photon-pion forward Compton scattering am- 
plitude. 7r, dashcd-linc; 7, wavy-line; 5, internal solid-line, dressed-quark propagator, Eq. 13.241 . 
The rilled circles represent the pion's Bethe-Salpeter amplitude, r„- in Eqs. 12.381 . l3.3UI - f3.33l . 
and the dressed-quark-photon vertex, r M in Eq. 13.351 . depending on which external line they 
begin/end. (Adapted from Ref. 1761 .'! 



forward Compton scattering that survive in the deep inelastic Bjorken limit: 

q 2 

q 2 — ► 00 , P ■ q — > —00 but x := fixed. (3-14) 

2P ■ q 

The upper diagram represents the renormalised matrix element 
,A 

T+(q,P)=tr J T_I\(fc r ;-P) 

x S(k t ) ieQT v (k t ,k) S(k) ieQT^k, k t ) S(k t ) r+r,(fc r ; P) S(k s ) , (3.15) 

where: t± — §(ri ± 1T2); S(£) — di&g[S a (£), Sd(i)], with S u = Sd — S, assuming 
isospin symmetry; and fcp = k—q — P/2, k t = k — q, k s = k — q—P. A new element in 
Eq. I|3.15|l is T tl (£i,£2), the dressed-quark-photon vertex, with Q = diag(2/3, —1/3) 
the quark-charge matrix. It can be obtained by solving the inhomogeneous vector 
BSEp-^i.e., Eq. I|2.9[) with 757^ — > 7^, or modelled, based on symmetry considera- 
tions, as we describe below. The matrix element represented by the lower diagram 
is the crossing partner of Eq. (|3.15() and is obvious by analogy. 

The hadronic tensor relevant to inclusive deep inelastic lepton-pion scattering 
can be obtained from the forward Compton process via the optical theorem: 

W^faP) = W+„(q;P) + W-„(q;P) = i-lm [T+ (g; P) + T" (q; P)] . (3.16) 
In the Bjorken limit one 

w+(q; p) = F+(*) V + F2 (x) 4^ , (3.17) 
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t^, v = <5 M „ - q^q v /q 2 , q] l = q l x + IxP^, and 

F+(x) = 2xF+(x) , F+ 2 (x) — > as x — > 1 . (3.18) 

Combining these results with their analogues for W~ , one recovers Bjorken scaling 
of the deep inelastic cross section, namely, the cross section depends only on x, 
and not on P ■ q and q 2 separately. Furthermore, the derivation shows that, in the 
Bjorken limit, x is truly the fraction of the pion's momentum carried by the struck 
quark, and one may therefore write 

F^{x) = F+{x) + F^{x) = ~[xu{x) + xu{x)] + ^[xd(x) + xd(x)} + (3.19) 

where u{x), u(x), etc., are the quark and antiquark distribution functions, and 
the ellipsis denotes contributions from the s— and c-quarks. (Heavier quarks are 
assumed not to contribute at all.) 

As we now explain, the calculation in Ref. [75j produces the valence-quark dis- 
tributions: 

qZ{x):=q*{x)-<?(x). (3.20) 

It is plain from Eqs. I|3.15|l . i|3.16|l that the hadronic tensor depends on the dressed- 
quark propagator, pion Bethe-Salpeter amplitude and dressed-quark-photon vertex, 
and although sea-quarks are implicitly contained in these elements, the "handbag" 
impulse approximation diagrams in Fig. l3.4l onlv admit a coupling of the photon to 
the propagator of the dressed-quark constituent. The quark's internal structure is 
not resolved. The calculation therefore yields the valence-quark distribution at a 
scale <7o that is characteristic of the resolution and £q = 1/qo is a length-scale that 
typifies the size of the valence quark. As with all calculations of this type hitherto, qo 
is an a priori undetermined parameter, although one anticipates 0.3 < qo < 1.0 GeV 
(0.7 > £q > 0.2 fm), with the lower bound set by the constituent-quark mass and 
the upper by the onset of the perturbative domain. A sea-quark distribution is 
generated via the evolution equations when the valence distribution is evolved to 
that g 2 -scale appropriate to a given experiment. An explicit sea-quark distribution 
at the scale qo can arise from non- impulse diagrams; e.g., contributions that one 
might identify with photon couplings to intermediate-state quark-meson-loops that 
can appear as a dressing of the quark propagator: 

7T+ = u d -> (uss) d = (K + s) d -> u d = tt+ , (3.21) 

with deep inelastic scattering in this instance taking place on the kaon: s — > 
K + sj, etc. Such intermediate states arise as vertex corrections in the quark-DSE 
and were neglected merely to simplify the first calculation of (x) , the improvement 
of which is a modern challenge. With these observations in mind, 

F+(x;q ) = lxu:(x;q ), F 2 " (x; q ) = | x %(x; q ) ■ (3.22) 

The calculations are required to yield 

dxul(x;q ) = 1 = / dx d£(x; q ) ; (3.23) 
Jo 
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viz., to ensure that the n + contains one, and only one, w-valence-quark and one 
d- valence-quark. 

3.3.2. Algebraic parametrisations 

To complete the calculation one must specify the elements in the integrand of Eq. 
(|3.15l) . This is where the algebraic parametrisations appear. The dressed- light- 
quark propagator is 

S( P ) = -i 1 - P av(p 2 )+<Js(p 2 )=[iJ-pA( P 2 )+B( P 2 )]- 1 , (3.24) 
a s (x) = 2mT{2(x + fh 2 ))+T{b 1 x)T{b 3 x) [b + b 2 T{ex)] , (3.25) 

a v (x) = —L^ [1 - T(2(x + to 2 ))] , (3.26) 
x + m z 

with J-(y) = (1 — c~ y )/y, x = p 2 /X 2 , fh = to/A, 0-5(2;) = \<Js(p 2 ) and <5y(ir) = 
A 2 oy(p 2 ). The mass-scale, A = 0.566 GeV, and dimensionless parameter values: 5 

m 6 h b 2 63 (3 27) 



0.00897 0.131 2.90 0.603 0.185 



were fixed in a least-squares fit to light-meson observablespSI and the dimensionless 
w-current-quark mass corresponds to 

m„(l GeV) = 5.1 MeV. (3.28) 

The pointwise form of the dressed-quark wave function renormalisation and mass 
function obtained with this simple parametrisation is qualitatively identical to that 
of the numerical DSE solutions depicted in Figs. 13.11 13.21 and hence is consistent 
with lattice-QCD simulations. However, it was proposed long before both, and pro- 
vided the first clear evidence that nonperturbative momentum-dependent dressing 
of parton propagators is fundamentally important in QCD and provides a key to 
understanding hadr on pro perties. The parametrisation represents the propagator 
as an entire functio n 1 1 and thereby exhibits confinement through the violation 
of reflection positivity,' 1 which means, loosely speaking, that the quark fragments 
before it can reach a detectorp^l and manifests DCSB with 

3 b 1 

- (OT)iGcV^ =^7Trf ln T2 = (0-221 GeV) 3 , (3.29) 
4tt bi b 3 A QCD 

which is calculated directly from Eqs. I|2.53|l . I|3.5|) after noting that Eqs. 13.251) . 
(|3~2"6l) yield Eq. l(2~4^|l with j m = 1. 



9 e = 10 4 in Eq. 13.251 acts only to decouple the large- and intermediate-p 2 domains. The study 
used Landau gaug e be cause it is a fixed point of the QCD renormalisation group and Z2 ~ 1, even 
nonperturbativelyED 

''This sufficient condition for confinement is discussed at length in Sec. 6.2 of Ref. 1211 . Ref. I8UI . 
Sec. 2.2 of Ref. |M] and Sec. 2.4 of Ref. [52]. 
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The general form of the pion's Bethe-Salpeter amplitude is given in Eq. 1(2.38(1 . 
The behaviour of the invariant functions therein is largely constrained by the axial- 
vector Ward-Takahashi identity, as we explained in connection with Eqs. (12.42(1 - 
((2.45(1 . These relations and numerical studied^ support the parametrisatiorP^I 



E n (k;P) = — B^k 2 ), 
where B n is obtained from Eqs. ((3.24|) - ((3.26|) . evaluated using fh = and 



bl = 0.19 



with the other parameters unchanged, and: 



F n (k;P) = EAk; P) / (no U); 
G„{k;P) = 2F 7r (k;P)/[k 2 + M 2 v ], 



(3.30) 

(3.31) 

(3.32) 
(3.33) 



Muv = IOAqcd; and H n (k;Q) = 0. The amplitude is canonically normalised 
consistent with the impulse approximation: 

<Pk 



(27T) 



tr C D 
+tr CD 



T^k;-P)^^T 7I (k;P)S(k^ 



dP^ 

r^(k;-P)S(k + )T 7r (k;P) 



dS(k. 
dP u 



(3.34) 



This fixes N v . The decay constant is subsequently obtained from Eq. 1(2.51(1 . 

The manner whereby an Abelian gauge boson couples to a dressed-fermion has 
been studied extensively and a range of qualitative constraints have been elucida- 
ted.ESl This research supports an Ansatz, 



8G 



iT^h,^) = iZA<AA)l* + (*i ■ (£i+£2)A A (e 2 1 J 2 2 ) + A B (el£ 2 2 )] , 

(3.35) 

wherein 



E F (£l,£ 2 2 ) = l[F(e 2 )+F(e 2 )}, A F (£lJ 2 2 ) 



(3.36) 



with F = A,B the scalar functions in Eq. ((3.241) . which preserves many of the 
constraints, in particular, the vector Ward-Takahashi identity. Furthermore, the 
Ansatz is expressed solely in terms of the dressed-quark propagator. 

Using the elements just described, one obtains the following calculated values 
for an illustrative range of pion observables (adapted from Ref. (84( ) : 





(GeV) 


/„• (GeV) rv(fm) 


Calc. 

Expt.™ 


0.139 
0.138 


0.090 0.56 
0.092 0.66 



(3.37) 



with the pion's charge radius, r T , calculated in impulse approximation!^ 
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x 



Fig. 3.5. Solid line: u v (x; go) calculated using the DSE model described in Sec. 3.3. The resolving 
scale go = 0.54 GeV = 1/(0.37 fm) is fixed as described in connection with Eq. 13.401 . Dashed line: 
u v (x;qo) calculated in the absence of the pseudovector components of the pion's Bethe-Salpeter 
amplitude; i.e., F = = G in Eq. 12.381 instead of Eqs. l3~32t and 13.331 . Dot-dashed line: 
distribution calculated with —* 0.1 m n , M = 0.36 GeV. (Adapted from Ref. 1761 . ) 

3.3.3. Calculated distribution function 

In the Bjorken limit, Eq. I|3.15|l reduces to a one-dimensional integral that depends 
parametrically on the valence-quark mass, M , and using the parametrisations de- 
scribed above that equation yields the valence-quark distribution function via Eqs. 

iprmji . (233, irrm The value 

M = 0.30 GeV (3.38) 

is fixed by normalisation, Eq. Ij3.23|l . and gives the distribution function depicted 
in Fig. 13.51 It vanishes at i = 1, in accordance with the kinematic constraint 
expressed in Eq. I|3.18|l . and corresponds to a finite value of Fi(x = 0), which 
signals the absence of sea-quark contributions. Unsurprisingly, since the pion is a 
light bound state of heavy constituents, the shape of the distribution is characteristic 
of a strongly bound system: cf. for a weakly bound systerr]SSIii t ,(a;) « S(x — |). 
The momentum-fraction carried by the valence-quarks at this resolving scale is 

( x qV — / dxx [u v (x; qo) + d v (x; qo)] = 0.71 , (3.39) 
Jo 

with the remainder, {x g Y — 0.29, carried by the gluons that bind the pion bound 
state, which are invisible to the electromagnetic probe.* As with all calculations of 
parton distributions hitherto, in Ref. |76j the resolving scale 

q = 0.54 GeV = 1/(0.37 fm) (3.40) 

'NB. The parametrised pionic parton distributions in Ref. 1891 yield a gluon momentum-fraction 
of (x g Y = 0.29 at go = 0.51 GeV. 
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Fig. 3.6. Dashed lineCSI xu v (x; <jo); short-dashed line: fit of Eqs. 13.421 and 13.461 . Solid line: 
the evolved distribution, xu v (x;q = 2GeV); dotted line: xu v (x;q = 4.05 GeV), evolved with a 
4-flavour value of Aqcd = 0.204 GeV; and dot-dashed line: the phenomenological fit of Ref. 1901 . 
The Drell-Yan data are from Ref. EZJ. (Adapted from Ref. 1751 .) 



was chosen so that when using the nonsinglet evolution equations (see, e.g., Ref. 

to evolve the distribution in Fig. 13. 51 up to q = 2 GeV, agreement was obtained 
between the first and second moments of the calculated distribution and those 
determined from a phenomenological fit to data, in this case the fit of Ref. HH, viz- 



LattP 



L g/2GcV 



t q/2GcV 



L q/2GcV 



0.24 0.098 0.049 

0.24 ±0.01 0.10 ± 0.01 0.058 ±0.004 

0.27 ±0.01 0.11 ± 0.03 0.048 ±0.020 



(3.41) 



The original and evolved distributions are depicted in Fig. 13.61 

A fit to xu v (x; q), acceptable for the estimation of moments, is obtained with 

iC^ <?) = A>,/3 x a (1 - xf, A a , = T(l + a + 0)/[T(a) T(l + (3)}, (3.42) 

whereupon the moments of u™ om (x; q) are given by 

OL-1 



n 



i + a + (3 



(3.43) 



The Drell-Yan data^I are described by Eq. (|3~^2l with q = 2 GeV values 

a DY = 0.57 ± 0.03 , /3 DY = 1.27 ± 0.04 , (3.44) 
while a global fit to Drell-Yan and prompt photon data yields the consistent resullP^I 
a fit = 0.64 ± 0.03 , /? fit = 1.15 ± 0.02 . (3.45) 
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Fig. 3.7. F uller comparison of theory with experim ent. Solid line: DSE resultpSl dot-dashed line: 
NJL modelE S^l l ong-dashed line: instanton modelfHU short-dashed line: light-front constituent 
quark modelEU All calculations have been evolved to q = 4.0 GeV using a 4-flavour value of 
Aqcd = 0.204 GeV. The Drell-Yan data arc from Rcf. EJ. (Adapted from Ref. 

The calculation of Ref. jZHl is described by the values 



q (GeV) 


0.57 


2.0 


4.05 




a 


1.34 


0.92 


0.84 


(3.46) 


P 


1.31 


1.80 


1.98 





A material feature of this DSE result is the value of j3 ~ 2 because although pertur- 
bative QCD cannot be used to obtain the pointwise dependence of the distribution 
function it does predidP^ the power-law dependence at x ~ 1: 

pQCD: <(a:) x x (1 - xf , (3.47) 

in agreement with the DSE result (corrections are logarithmic). However, as is 
apparent in Fig. 13.61 and emphasised by Eqs. Q3.44p . 13.45|l . this prediction disagrees 
with extant experimental data.-' That is very disturbing because a verification of 
the experimental result would present a profound threat to QCD, even challenging 
the assumed vector-exchange nature of the force underlying the strong interaction. 

In Fig. 13.71 we illustrate that the only extant calculation which agrees with the 
distribution inferred from irN Drell-Yan data is that performed with a particu- 
lar regularisation of the Nambu-Jona-Lasinio model!^ That calculation yields a 
distribution 

uf(x; q^ JL = 0.35 GeV) = 6{x) 0(1 - x) (3.48) 

J This is in spite of the fact that the first four moments agree: Eq. 13.411 . Plainly, the low moments 
contain little information about the distribution on the valence-quark domain. Lattice-simulations 
are currently limited to the low moments. 
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Fig. 3.8. Deep inelastic scattering from the proton's 7r-cloud (Sullivan process) could provide a 
means of measuring u^(x) at JLab. (Adapted from Rcf. 1961 .) 

which corresponds to the valence-quark carrying each and every fraction of the 
pion's momentum with equal probability. The result is an artefact tied to this 
model's representation of the pion bound state by a momentum- independent Bethe- 
Salpeter amplitudef^ESI i.e^ representing the pion as a point-particle, which is 
a necessary consequence of the model's momentum-independent interaction. It is 
clearly the hardest distribution that is physically possible. 

This observation underscores the serious nature of the discrepancy described 
above in conn ection with Eq. (|3.47l) . and in highlighting that disagreement the 
DSE studjC^has catalysed interest in u*(x) and proposals for its remeasurement. 
One proposal that could use existing facilities would employ the process depicted 
in Fig. 13.81 at JLabP^ This process could also be used efficaciously at a future 
electron-proton collider to accurately probe u%(x) on the valence-quark domainPSI 

3.3.4. Summary 

This application illustrates the power of using algebraic parametrisations of key DSE 
elements: one can proceed rapidly to an insightful analysis of hadronic phenomena 
and arrive at robust conclusions. Now, however, this type of analysis is supported, 
improved and in some cases superseded by a direct application of the one-parameter 
rainbow-ladder model specified by Eq. (|3.1() . 

4. Ab Initio Calculation of Meson Properties 

In this section we review the systematic application of the renormalisation-group- 
improved rainbow-ladder kernel defined by Eq. (|3.1|) to a wide range of meson 
observables. We emphasise at the outset that each result is a prediction, in the 
sense that the model's mass-scale is fixed, Eq. (13.311 . and every element in each 
calculation is completely determined by, and calculated from, that kernel. 

4.1. Elastic electromagnetic form factors 

We first consider processes that involve only three external particles, at least one of 
which is a hadron. The simplest of these are the pion and kaon electromagnetic form 
factors. The matrix element describing the coupling of a photon with momentum q 
to a pseudoscalar meson composed of quark a with electric charge e a and antiquark 
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b, electric charge eg, can be written as the sum of two terms 

A? (P, q) = e a A? a (P, q) + e 5 Af\P, q) , (4.1) 

which, respectively, describe the coupling of the photon to the quark and the anti- 
quark. In the Breit frame, where the incoming meson has total momentum P — q/2, 

A ab (P,q) = 2P v F{q 2 ), (4.2) 

where F{q 2 ) is the meson's elastic electromagnetic form factor, (r 2 = — 6F'(0) is 
the square of the meson's charge radius.) Each term on the r.h.s. of Eq. 1)4. lfl has 
this property and hence one can also write 

Af(P,q) - 2P„[e a F a (q 2 ) + e- b F b (q 2 )}, (4.3) 

and thereby explicate the contribution of each quark to the total elastic form factor. 

The matrix element describing the scatter ing process is most easily calculated 
using the renormalised impulse approximation, in which, e.g., 

,A 

Af\P,q) =iV c tr D / S a (£)T ab (r + ;P^)S b (£ + )iTl(£ + ,£^S b (£„)T ba (r_;-P + ), 



(4.4) 

where: £ = k + P/2, t± = k - P/2 ± q/2, r± = k ± q/A, P ± = P ± q/2; and, 
similar to what we saw with Eq. 1)3. 15)1 . S a ' b is a dressed-quark propagator, T ab is 
the meson's Bethe-Salpeter amplitude, and r^ ,b is a dressed-quark-photon vertex. 

In this section the dressed-quark propagators are obtained by solving the mass- 
dependent rainbow gap equations, Eq. (|2.13(l . and the Bethe-Salpeter amplitudes 
come from Eqs. 1)2.12(1 . 1|2.37J) . Furthermore, the dressed-quark-photon vertex is 
calculated by solving a renormalised inhomogcncous BSE, viz. 

r£(fc + ,fc_) = Z 27M + / K(k,£;q)S a (e + )T;(£ + ,£„)S a (£-), (4.5) 
J i 

wherein k± — k ± q/2, £± — £±q/2, with the kernel again given by Eq. I|2.12|l . At 
this point it is important to remember that because the rainbow-ladder truncation 
is the first term in a systematic DSE truncation scheme [Sec. 2.2] the vertex thus 
obtained satisfies the Ward-Takahashi identity 

i^rjC^,*-) = s//\k + ) - s-\k.) . (4.6) 

(NB. This is only true if a Poincare covariant regularisation scheme is employed.) 
4.1.1. Current conservation 

It is electromagnetic current conservation that reduces to one the number of form 
factors on the r.h.s. of Eq. 14. 2|) and it also requires F a (0) = 1 = F b (0). As with the 
realisation of Goldstone's theorem, these consequences of symmetry are exhibited 
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Fig. 4.1. Top row: The three additional contributions to the normalisation of the Bethe-Salpeter 
amplitude, Eq. 13.341 . generated when the second term in Eq. 12.141 is included in the gap equa- 
tion's kernel. (Without loss of generality, for this illustration the integration variable is chosen 
such that the total momentum, P, flows only through the antiquark line, upon which the slash rep- 
resents dS/dPfx.) Bottom row: The corrections to the impulse approximation, Eq. 14.41 . necessary 
to maintain current conservation. (Adapted from Ref. 1371 .) 



by the electromagnetic matrix elements if, and only if, there is an intimate relation 
between each of the elements in the calculation. The impulse approximation, calcu- 
lated with propagators and vertices obtained using the rainbow-ladder truncation, 
preserves these constraints without fine-tuningp"2l i.e., independent of the detailed 
form of the interaction, Eq. . 

Suppose now that one adds the second term of Eq. (|2.14|) to the kernel of the 
gap equation. Following the systematic procedure reviewed in Sec. 2.2, that term 
generates three additional contributions to the kernel of the BSE and also intro- 
duces a dependence on the bound state's total momentum. That in turn gives rise 
to a modification of the canonical normalisation condition for the Bethe-Salpeter 
amplitude; viz., ° the integrals represented by the three diagrams in the top row 
of Fig. 14.11 must be added to the r.h.s. of Eq. I|3.34fl . The dressed-quark-photon 
vertex calculated from Eq. I|4.5II with the modified kernel automatically satisfies the 
Ward-Takahashi identity, Eq. Ij4.6|l , and current conservation is guaranteed without 
fine-tuning provided one augments the impulse approximation by the three dia- 
grams depicted in the bottom row of Fig.ETTl 37 That these diagrams are necessary 
and sufficient is obvious once one realises that, according to the differential Ward 
identity: iT^(k,k) — dS(k)/dk fJ- , the derivative of a quark line is equivalent to 
the insertion of a zero momentum photon. The generalisation of this procedure is 
straightforward and it is therefore apparent that for any given gap equation kernel 
it is systematically possible to construct the current-conserving matrix element that 
describes the coupling of the photon to a composite meson. 
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4.1.2. Pion and kaon electromagnetic form factors 

To be concrete we review a calculation of the ir and K form factorsl 



FAq 2 ) = \k{i 2 ) + \f^), 
F K «(Q 2 ) = 4^o(? 2 ) + ^ (<z 2 ). 



(4.7) 



(4.8) 



(4.9) 



So far as the strong interaction is concerned u and d quarks are identical in the 
isospin symmetric limit and hence there are only three independent form factors, 



F^(q 2 ) = F*{q 2 ), F£ + {q 2 ) = Ff <a (q 2 ), and F s K+ (q 2 ) = F s K0 (q 2 ), which in impulse 



approximation are given by Eqs. (|4.3|1 . Ij4.4|l . 

These form factors were calculated in Ref. )37| as an ab initio and parameter- 
free application of the model defined by Eqs. (|3.3[) . Therein the dressed-quark 
propagators were calculated and used in constructing the kernel of the BSE; the 
BSE was solved to determine the mesons' bound state amplitudes completely, i.e., 
all the amplitudes in Eq. (|2.38|) were calculated and shown to play an important 
role; and then the elements were combined to yield a manifestly Poincare covariant 
prediction of meson electromagnetic form factors. 

The result for the pion form factor is depic ted in Fig. l4.2l fc wherein it is compared 
with the most recent experimental dataP^I The pion's calculated charge radius is 
compared with experiment in Table 14.11 The excellent agreement is misleading 
because it-it final-state interactions add < 15% to the impulse approximation result 
for tvI ^^ I (NB. This rescattering contribution to the form factor vanishes with 
increasing q 2 .) However, that only makes the result more plausible by leading to a 
disagreement with experiment commensurate with that one would anticipate based 
on the discussion in Sec. 2.2. 

It is evident in Table 14.11 th at th e calculated charged kaon form factor, Eq. 
(|4.8|1 . agrees with availabielata,™ which unfortunately is not a stringent test 
of theory because it covers only a small low-g 2 domain and has large errors. The 
g 2 -dependence of the form factor again deviates from a simple monopole on q 2 > 



The neutral kaon form factor, Eq. (|4.9I) . as the difference between two terms, 
reacts most to a model's details. Of course, F K o(q 2 = 0) = as a direct result 
of current conservation and Ref. |37j is guaranteed to reproduce that. However, 
the evolution away from q 2 — is sensitive to differences between the dressed d- 
and s-quark mass functions. (This is obvious because if they were the same then 
F K o = 0, just as F n o =0.) The calculation yields r 2 K0 < 0, Table ETT1 which is 
readily understood: the heavier and positively charged s-quark in the K° is more 

fc The nature and meaning of vector dominance is discussed in Sec. 2.3.1 of Ref. Iritil , Sec. 2.3 of 
Ref. 1811 and Sec. 4.3 below: the low-q 2 behaviour of the pion form factor is necessarily dominated 
by the lowest mass resonance in the J PC = 1 channel. Any realistic calculation will predict 
that and also a deviation from dominance by the p-meson pole alone as spacelike-g 2 increases. 



2GeV 2 . 
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Fig. 4.2. Impulse approximation DSE prediction for q 2 F 7T (q 2 ) obtained in an ab initio, parameter- 
free application of the renormalisation-group-improvcd rainbow-ladder truncation, Eqs. 13.11 . 
<3"3t . The data are from Refs. kWI^HMI (Adapted from Ref. ITo*2l 1 



often to be found at the meson's core so that the lighter <i-quark provides an excess 
of negative charge at the surface. The resul t is therefore qualitatively reliable. 
However, as the squared-charge radius is smallP^ K°-K° final state interactions, 
neglected in Ref. [37] . can have more of an impact on this observable. These effects 
disappear with increasing q 2 but data a re d ifficult to obtain for q 2 ^ 0. 

It is a model independent predictiorl^l of the DSE framework reviewed herein 
that the elastic electromagnetic meson form factors display 

q 2 F(q 2 ) = constant, q 2 > Aq CD , (4-10) 

with calculable ln d Q 2 /Aq CD corrections, where d is an anomalous dimension.' This 
agrees with earlier perturbative QCD analyses! ^7 | lU8 l However, to obtain this re- 
sult in covariant gauges it is crucial to retain the pseudovector components of the 
Bethe-Salpeter amplitude in Eq. l|2.H8fl : F n , G v . (NB. The quark- level Goldberger- 
Treiman relations, Eqs. 12.43JI . I)2.45[l . prove them to be nonzero.) Without these 
amplitudes}^ q 2 F(q 2 ) oc l/q 2 . Similar statements are true of the role played by 
nonleading components in the Bethe-Salpeter amplitudes of vector mesons. The 
calculation of Ref. |84j , which uses the propagator and vertex parametrisations de- 
scribed in Sec. 3.3, suggests that the perturbative behaviour of Eq. H4.10[l is unam- 
biguously evident for q 2 > 15 GeV 2 . Owing to challenges in the numerical analysis, 
the ab initio calculations of Ref. |37j cannot yet make a prediction for the o nset of 
the perturbative domain but progress in remedying that is being made^^ 

1 Corrections to the rainbow-ladder-impulse approximation, such as those depicted in Fig. 14.21 
contribute to the anomalous dimension but do not modify the power law dependence. 
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Table 4.1. Comparison of calculated ^7 l ? 7 lsquared -charge-radii (in fm 2 ) with datalHZnnSESlThe 
7*7r°7 interaction radius is also includedl I ' (Adapted from Ref. 1371 . ) 





r 2 

' 7T 


r 2 


r 2 


r 2 

7T77 


Calc. 


0.45 


0.38 


-0.086 


0.41 


Expt. 


0.44 ±0.01 


0.34 ±0.05 


-0.054 ±0.026 


0.42 ±0.04 


Rel. -Error 


0.023 


0.118 


0.59 


0.024 



4.2. K13 decays 

Related to the elastic form factors are the semileptonic transition form factors de- 
scribing K+ -> n £isi [K+], K° -► -n-lvt [K° a ] and vr+ -► 7r e^ e [7r e3 ]. They 
proceed via the flavour-changing vector piece of the V — A electroweak interaction, 
in particular j and . The axial- vector component does not contribute because 
the two mesons involved have the same parity. Neither nor j'^J" is conserved 
and the symmetry breaking term measures the dressed-quark mass difference [see 
Eq. I|4. 131)1. The se processes can therefore be employed to probe flavour symmetry 
violation j ^W^ ' ■ ' I and this more effectively than the neutral kaon form factor. 

The matrix element for the K® 3 transition is characterised by two form factors: 

Ai su ( P K,q,-p,) = (n-{p,)\- S ^u\K\p K )) = 2P p /fV) +9p/^V) , (4.11) 

where p^ = P — q/2, px = (P + q/2) and q is the VT-boson's momentum, with 
t = —q 2 . In the isospin symmetric case considered in Refs. |70lll()j . m u = and 
j£ u is conserved so that /£(*) = -F w (t), fl = 0. In addition, /f° = f£ + and 
hence all new information is contained in the two form factors of Eq. (|4. lip . Their 
calculation in impulse approximation involves only one element not already used 
in Sec. 4.1.2, namely, the dressed-suVF-vertex, which replaces the dressed-quark- 
photon vertex. That piece which contributes to decays is obtained from 

r»(fc + ,fc_)=Z 27|i + I K{k,b,q) S s {l + )T™(l + ,l_) S u (t-) , (4.12) 
Je 

where again, for a consistent truncation, the renormalisation-group- improved ladder 
kernel must be used. The vertex thus calculated satisfies 

i^r»(fc+,fc_) - S-\k + ) - S-\k^) - (m s - m„)rf (k+,k-) , (4.13) 

where r| u is a flavour-dependent scalar vertex analogous to the pseudoscalar vertex 
in Eq. (|2.48l) . The parameter-free predictiorlHSI for f? is depicted in Fig. 14.31 

On the physical domai n for decays: m 2 < t < {rax — m^) 2 = 0.13 GeV 2 , 
the calculated form factorJli^l are well approximated by a straight line: 



A 



mif'(mi)/f(0) , (4.14) 



supporting an assumption common in the analysis of experimental data 



EH 
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Fig. 4.3. DSE results for the K l3 form factor f + (t = -q 2 ). Experimental data from CPLEArIM 
(Adapted from Ref. IHol .1 



The last term in Eq. I|4.13|l indicates the manner in which decays are sensi- 
tive to the c urre nt-quark-mass difference and its enhancement through nonpertur- 
bative effectsPHl Experimentally this is best explored via the scalar form factor 

fo(q 2 ) = U(q 2 ) / 2 f-{q 2 ) , (4.15) 
m K — 

which measures the divergence q fl A'lf u (pK, q, — p^)- Current algebra predicts the 
value of fg- at the Callan-Treiman point jH^I t = m 2 K — m 2 =: A f$ (A) = 
—fa /fir = —1-23, while a systematic analysis of corrections yields! I f^(A) — 
— 1.18. The Callan-Treiman point is not experimentally accessible but the robust 
nature of its derivation makes the value of /o(A) a tight constraint on any theoret- 
ical framework. 

In Table POl we display a representative sample of the results calculated in Ref. 
(1 lUj . We no te that |/+(0)| « 1 and that is consistent with the Ademollo-Gatto 
theorem) -*- which states that flavour symmetry breaking effects are suppressed 
at t = 0. Furthermore, it was observed in Ref. |115j and confirmed in Ref. |7U|. 
that /-(0) is a measure of the s:u constituent-quark mass ratio. Applying this 
notion to Ref. in which the calculated value of that ratio is ~ 1.25, one 

estimates /-(0) « —0.15, from Fig. 1 of Ref. |115| . The calculated value is actually 
/_(0) — —0.10. These observations again emphasise the power of a systematic 
symmetry preserving DSE truncation. 

4.3. Form factors in the timelike region 

The electromagnetic pion form factor, F 7T (Q 2 ), exhibits a peak at q 2 = —t = —m 2 
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Table 4.2. Calculated^ A";3 observables, compared with experimental data" 3 '-' and chiral pertur- 
bation theory^Sl The partial widths are calculated in the usual way. (Adapted from Ref. ITTol .1 





DSE 


Expt. 


ChPT 


-/+(o) 






0.96 




0.98 


-/o(0) 






0.10 




0.16 


-/o(A) 






1.18 




1.18 


Af 






0.027 


0.0300 ± 0.0020 


0.031 


AT 






0.027 


0.0282 ± 0.0027 


0.031 


A^ 






0.018 


0.030 ± 0.005 


0.017 








0.018 


0.013 ± 0.005 


0.017 


r(^ 3 )xio- 


8 cV 


0.49 


0.494 ± 0.005 




r(^ 3 )> 


<io- 


-8 oV 


0.32 


0.346 ± 0.004 




rK)> 


10" 


8 cV 


0.24 


0.259 ± 0.003 




r(^+ 3 )> 


(10" 


-8 oV 


0.16 


0.174 ± 0.003 





that is associated with e + e — > p — > it + it . This is a general feature of the elec- 
troweak form factors of pseudoscalar mesons, all of which can be expressed as 

F P {t) = gvpp — ; = (4.16) 

m v — imvi v — t gv 

in the neighbourhood of the relevant flavour channel's lowest-mass vector meson 
resonance, assuming it is narrow. In Eq. (|4.16|) . gvpp is the coupling constant 
modulating the V — * PP decay, Ty is the total width of the vector meson and 
ray I gv characterises the strength of the vector- meson-photon coupling. 

This behaviour arises naturally in applications of the rainbow-ladder kernel. In 
the case of the pion and kaon form factors the poles appear as a straightforward 
consequence of solving the inhomogeneous BSE, Eq. (|4.5|l : a result we have already 
illustrated in connection with the axial-vector vertex, Eq. 112.4011 . In the neighbour- 
hood of the pole the solutions of Eq. (|4.5(l assume the formal 

rj(p+,P-) w fvm v ^— Tf{p;q) , (4.17) 

where is the vector meson's Bethe-Salpeter amplitude and fy was introduced 
in Eq. (|3.7|) . m The absence of a width in Eq. (|4.17|) is a defect of the rainbow-ladder 
truncation, which can be remedied without overcounting by including diagrams in 
the kernel that represent the vector meson's coupling to the PP intermediate state. 
Alternatively, if th e width:mass ra tio is small, that can be done via bound state 
perturbation theory! 74 | 1 ie> | ll 7 1 H8 | 

The timelike behaviour of the form factors describing electroweak 7r and K tran- 
sitions is depicted in Fig. 14.41 It is apparent that is singular at t = (0.74 GeV) 2 
and Ft. + at t = (1.05 GeV) 2 ; namely, at the masses of the p- and 0- mesons (see 

m fv and gv are algebraically related; e.g., for the p°, f p m p = \pivn? jg p \ and f^m^ = 3m'i/g^. 
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Table . Since every other element has already been calculated, the coupling 
constants g^uu, 94>kk can be determined by fitting Eq. i|4.16[l in the neighbourhood 
of the pole! I Alternatively, one can directly calculate the coupling constant s us- 
ing impulse approximation three-hadron triangle diagrams. The results coincid 
and are listed in Table l4~3l The agreement with experiment at the 9% level is good, 
given the rainbow-ladder truncation's established accuracy. 

Figure 14.41 also reveals poles in the electroweak transition form factors. 
The transverse form factor, f+(t), exhibits a pole at t = (0.94 GeV) 2 ; i.e., at the 
mass of the K* meson (see Table 13.211 , as expected, and it is this singularity that 
facilitates the successful description illustrated in Table !4.2p Sl The scalar form 
factor, fo(t), possesses a pole at t = (0.89 GeV) 2 . As the discussion makes clear, 
this is the model's prediction for the mass of the lightest 0+ s meson. The model 
also yields a 0+_ d ^ meson with mass 0.67 GeV.1^41 However, as observed in Sec. 
2.2, the rainbow-ladder truncation is known to be unreliable in the scalar channel 
and, until a quantitative analysis of the corrections is complete, it is impossible 
to associate these dressed-quark-antiquark bound states with any mes ons in the 
hadron spectrum or with the scalar states in model studies 
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Table 4.3. Calculated coupling constants for the two-pseudoscalar decays of vector mesons. The 
rms rel. -error is 9%. (Adapted from Ref. 11201 . ) 



Calc. 


5.2 


4.3 


4.3 


4.1 


Expt. 


5.99 ±0.02 


4.48 ±0.04 


4.58 ±0.05 


4.59 ±0.05 


Rel. -Error 


0.13 


0.04 


0.07 


0.11 
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Table 4.4. C alculat ed I -^ fl jjj co upling constants for radiative vector meson decays and associated 
widths compared with datal^l The calculations yield the ratios gyp 7 /my , which are tabulated: 
the rms rel. -error is 12%. Errors in the calculated values of meson masses propagate into the values 
of the widths [see Eq. <4"T3l l (Adapted from Ref. 



^ (GeV- 1 ) 


p — > 7T ± 7 


LU — > 7T7 


K*± -> ^±7 




Calc. 


0.69 


2.07 


0.99 


1.19 


Expt. 


0.73 ±0.05 


2.35 ±0.07 


0.84 ±0.05 


1.27 ±0.07 


Rcl. -Error 


0.055 


0.119 


0.18 


0.063 


r y ^ P7 (keV) 


Calc. 


53 


479 


90 


130 


Expt. 


67 ±7.5 


734 ± 35 


50.3 ±4.6 


117 ± 10 



4.4. Vector meson transition form factors 

Vector meson transition form factors, such as 7*7r — > p and 7*7r — ► lu, are important 
in the description of hadron photo- and electro-production reactions, and their 
analysis can help in developing a deeper understanding of the relation between 
QCD and efficacious meson-exchange models of the N-N interaction! ^^ A large 
variety of these form factors have been calculated in impulse approximation j E^ I ^ 1 
which for the p — > ttj transition is 

A£7(P;g) = ^tx Ti J A S(q 2 )r^(r + .-,-P-q)S(q 1 )T^r +0 -,P)S(k)irZ(k,q 2 ) 

= e^e^ al3 P a q F p ^(q 2 ), (4.18) 

TYlp 

where P is the p-meson's momentum, q is that of the photon and q\ = k + P, 
q2 = k — q, r a p = k + aP/2 + (3q/2. The coupling gp^-y is defined such that 
Fpn-yiq 2 — 0) = 1 and then 

2 m * 



r^=-^«^m p ^l-^] , (-1.1!)) 



where a cm is QED's fine structure constant. For m u = m^, T„±^ n ±„ = T p , 



0_> T 0^ 



and r w ^7r 7 = 9T p ^, n ^. The analys is of K* — > K^j is also straightforward and in 
the limit of S'[/(3)-fiavour symmetry! ^ 2 ! T K *±^ K ± 1 = r p ±^ 7r ± 7 and Yk*°^k°i = 
4Tjc*±_ ^±~, so one anti cipates r 'k*°-*k i > ^K k± ^K ± j m reality. 

The calculated result J 102 ! 1 ^ are summarised in Table 14.41 The agreement 
is consistent with the established accuracy of the rainbow-ladder truncation. The 
transition form factors are d epicted in Fig. l4.5l and. on the physical domain, F p7T1 (q 2 ) 
agrees with extant dat eP^ to a degree consistent with the absence of a width for the 
/3-meson in rainbow-ladder truncation. An algebraic analysis of the behaviour of the 
form factors at large spacelike-g 2 has not been completed but they do appear to fall 
faster than l/q 2 (see the w-dominance curve in the figure). This is consistent with 
analyses in perturbative QCD, which indicate that such processes are suppressed 
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Fig. 4.5. Left panel: DSE result for F UJW1 (q' 2 ) = F p7r -/(q 2 ) co mpare d with experimental datali^l 
The dot-dashed line is an earlier, rudimentary DSE calculations!" proceeding from parametrisa- 
tions of the propagators and vertices after the fashion described in Sec. 3.3. Right panel: DSE 
result for Fk *K-y{q 2 ), charged (dashed line) and neutral (dotted line). The solid line is the ratio 
F Kt+K+ (q 2 )/F K {q 2 ). (Adapted from Refs. 110511581 .1 



because they do not conserve hadron helicity.^3 The charged K*Kj form factor 
falls very rapidly because of cancellations between the contributing j-u- and j-s- 
quark diagrams. The analogous contributions for K*° interfere constructively. 

4.5. 7T° — > 77 

The decay 7r — > 77 and associated 77* — » n° transition form factor are closely 

related to the processes discussed in the last subsection and furthermore they have 

long been recognised to p osses s a number of unique features that are especially 

1 OS 

important for testing QCD ! UQ | 

In impulse approximation the transition form factor is obtained from the vertex 

= a em ^^N c truj S(£ 1 )T*(£;-P)S(£ 2 )iT^£ 2 ,£ 12 )S(£ 12 )ir„(£ 12 ,£ 1 ) 



= 2^ em _ J* 11 e^pa- kipk 2c r Fny-y (h ] fc 2 ) , (4.20) 



where l\ = £ - ki, 4 = ^ + ki, f- = (h + 4j)/2, £i 2 = £ - h + k 2 , and the coupling 
constant g W77 is defined such that F nj ^(ki = = k?) = 1. For k\ = = k 2 , the 
vertex describes 7r° decay and yields the width 
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2 /[8tt 2 / 2 ]). The data are from the CELLcEHO and CLEdSOl 
collaborations (Adapted from Rcf. 111)51 .) 



a monopole with mass scale m 2 
asymptotic form 1/(1 



It is a textbook example that g 7ril = and tt — ► 77 is forbidden in the absence 
of the Abelian anomaly; i.e., if the axial-vector quark current J| M = g| ^?l5l^<l 
is conserved in the chiral limit. However, this current is anomalous, as may be 
demonstrated in many ways, and one arrives instead at the model-independent 
chiral limit result: = ±, wherewith Eq. (|4.21(l predicts T^o^ = 7.73 eV cf. the 
experimental value!^ 7.84 ± 0.56, which corresponds to = 0.504 ± 0.018. 

It is a key feature of the DSEs that, with a systematic and nonperturbative 
truncation scheme, all consequences of the Wess-Zumino term and Abelian anomaly 
are obtained exactly, without fine tuning. A true representation of DCSB and the 
preservation of Ward-Takahashi identities is crucial in ach ieving this! ^^^^^ 
Hence, one obtains algebraically g^ l (m 1T = 0) = A anc J105| 

5*77 Kr) = 0.502. (4.22) 

The 7* 7r — > 7 transition form factor, which was measured by the CELLcPSI 
and CLEcPSni 

collaborations, is defined as 

F^^iq 2 ) = F^(kl = q 2 , k 2 = 0); (4.23) 

i.e., as the form factor in Eq. (|4.20JI evaluated with one photon on-shell, and the 
impulse approximation result is plotted in Fig. I4.6I The calculated interac tion 
radius r 2 7J := — 6i^, 7r7 (0) = 0.41 fm 2 agrees with the experimental estimateP^H 
0.42 ± 0.04. It is evident that the constraints of DCSB are very tight. 

The behaviour of the form factor at large spacelike-g 2 has been analysed in 
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perturbative QCD,"with the leading order result PHEHHIIIHS 



i^V) = J 0)-f- + O\-^-,-^y (4.24) 

J{w) = U dx i- W H2 X -ir (4 - 25) 

where w = (kf — + ^1) = 1 f° r l* n ~ * 7' an d ^(^j m( z) is the pion's light- 

cone quark distribution amplitude, J Q dx (/^(x; lnq) = 1. In this context accurate 
measurements of the transition form factor can be interpreted as constraining the 
x-dependence of 4> v (x;hiq) at the experimental scale, lng. At very large q 2 ; viz., 
ln[<7 2 /AQ CD ] ^> 1, ^(x; oo) = 6x(l — x) and consequently J(l) = 2. Existing data 
at ln[g/A QCD ] ~ 2 favouJ^OI 

J(l) w 1.5, a value which requires either <^> 7r (x;~ 2) 
to be far from asymptotic (e.g, <f>„(x;~ 2) = 630x 4 (l — x) 4 gives J(u>) = 1.5) or 
material corrections to the leading order result EZH 

The expression in Eq. (|4.20() has also been used to analyse th e large spac elike-q 2 
behaviour of the 7*7r — > 7 transition form factor, with the resullP2C231I25l 

4 4tt 2 f 2 

iW<? 2 )= o" T^- (4.26) 
o q 

This was obtained by assuming that for k 2 = and k\ — q 2 — > 00, 2fci ■ &2 = — Q 2 
because P = k\ + &2 and the pion is on-shell, and 

S(4 2 ) = i- -4 h ■ {ki - fe) , r M (i 2 , 4 2 ) = 7^ = • (4-27) 

The effect of Eq. (14.27(1 is to treat the photons symmetrically in the integrand of Eq. 
(|4.2U|) . In reality that might effect the situation k\ = q 2 = k 2 , which corresponds to 
W = 0. Since J(w = 0) = 4/3, this interpretation would reconcile Eqs. I|4.24|l and 
(|4.26|> . In essence, this expresses the relevant observation in Ref. |136| and suggests 
that Eq. I|4.20|l should be reanalysed with a more conscious focus on the asymmetric 
7*7r — > 7 transition. 

The calculations reported in Ref. |105j are not affected by these considerations. 
They simply provide the prediction of the rainbow-ladder DSE model for _F' 7 » 7r7 (g 2 ) 
on the domain covered by extant experiments, which do not extend into the truly 
asymptotic domain, ln[g 2 /Aq CD ] 3> 1. The result is plotted in Fig. 14.61 and on the 
scale of the figure is well approximated by a simple p-meson monopole. It agrees 
with existing data, within errors, and interpreted in terms of Eq. (|4.24() . corresponds 
to J(w = 1;1hq/A Q cd ~ 2) « 1.7. 

4.6. Scattering processes 

Consider now the more complex situation of four external lines. We have repeatedly 
emphasised that the DSEs provide for a resolution of the dichotomy of the pion as 

n NB. The anomalous dimension is zero; i.e., at leading order, in contrast to the elastic electromag- 
netic pion form factor, F 1 * 1T1 does not depend on the strong running coupling, a B (q 2 ). 
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Fig. 4.7. Sum of terms required to describe tt-tt scattering when the rainbow-ladder truncation 
is used to calculate the dressed-quark propagator and pion Bethe-Salpeter amplitude. In every 
case, the interaction line is specified by Eq. 13.11 . The first diagram on the top line, Do(s,t,u), 
is expressed in Eq. 14.281 : the remaining sum on the top line is denoted by Dq- S (s,t,u); and the 
sum on the bottom line by Dq- t (s,t,u). (Adapted from Ref. liitil .) 



both a Goldstone mode and a bound state of massive dressed-quarks, and the role 
played in this by the systematic, nonperturbative truncation scheme reviewed in 
Sec. 2. As our exemplar we therefore choose tt-tt scattering. It will immediately be 
appreciated that this process has long been of interest and, i ndee d, that its partic- 
ular features and the effective Lagrangian that describe theirP^I are a keystone of 
chiral perturbation theory. 

Based on the analyses reviewed in the preceding subsections one would naively 
suppose that the dominant contribution to tt-tt scattering arises from the first dia- 
gram in Fig. 14.71 which represents the renormalised expression 

A>(s, t, u) = tr / r«(£; -p 4 ) 

J i 

x r^ ++0 ; -ps) S(e +++ ) T^i +0+ ; P2 ) S(i+-+) r„(£ _ +;Pl ) , (4.28) 

with s = -(pi + p 2 ) 2 , t = -{pi - p 3 ) 2 , u = -(pi - p^ 2 , pi +p 2 = Pa + Pi and 
(-api — £ + fPi + §P2 + %P3- However, it follows immediately from Eqs. I|2.38|) . 

jsrsn tnat in tne cniral limit 

0, 0) = tr isosp ^VW] ^ I P^ff^p (4-29) 

and hence D o (0, 0,0) ^ when chiral symmetry is dynamically broken. It is clear 
now that this contribution alone violates the current algebra results for low-energy 
it-it scattering; viz., at threshold the scattering amplitude must vanish as m 2 /f 2 . 

The observation is not new. In fact, it has long been known that in the chi- 
ral limit D o (0,0, 0) is cancelled by contributions that may be described as scalar- 
isoscalar two-pion correlations. That mechanism is readily realised; e.g ., by analysi ng 
the auxiliary field effective action in four-fermion interaction theoriesl^SIESISSl It 
can also be achieved directly using a systematic truncation of the DSEs, by which 
manner it is firmly placed in context with the material reviewed herein. 
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For the latter, if one employs a dressed-quark propagator and Bethe-Salpeter 
amplitude obtained in the rainbow-ladder truncation, then the combination of di- 
agrams depicted in Fig. 14.71 is gu arantee d to reproduce the current-algebra results 
for near-threshold it-it scattering PSC2SI To illustrate, we note that the sum in the 
figure can be expressed 

D(s,t,u) = D (s,t,u) + D Ts (s,t,u) + D Tt (s,t,u) (4.30) 

and, in the chiral limit, it is readily established" that if any one of the external pion 
momenta vanish then 

D Ts (0,0,0) = -AA>(0,0,0) = D % (0,0,0) (4.31) 

and thus 25(0,0,0) = 0. Extending this, one finds at leading order in powers of 
mass andmomentM 

1 

from which follows the isospin-zero, -one and -two scattering amplitudes: 

2s — ml 



D(s,t,u) = —^(s + t-u), (4.32) 



Ao(s,t, u) — 3 D(s,t,u) + 3 D(s,u,t) — D(u,t, s) = 

t — u 

A 1 (s,t,u)= 2D(s,t,u)-2D(s,u,t) = -r-pr , (4.33) 



A 2 (s,t,u)= 2D(u,t,s) 



— 'i 

K - 



viz., precisely the current algebra results for low-energy it-it scattering, independent 
of the detailed form of the interaction. 

In Fig. 14. 81 we plot the mass-dependence of the isospin-zero and -two scattering 
amplitudes at threshold. At the physical value of the light current-quark-mass, 
95% of D (Am 2 , 0, 0) is cancelled by the remaining terms in Fig. 14.71 Furthermore, 
the direct calculation yields Aq = 8.6, A^ = —2.3 cf. the current-algebra results: 
A^y = 7toJ/(2/2) = 7.8, A 2 W = -m\jfl = -2.2; i.e., the complete result exhibits 
a 10% increase over the leading order value in the isospin-zero channel, which is an 
harbinger of the resonant behaviour of tt-tt scattering in this channel. 

For a detailed comparison with experiment, as with meson charge radiifi^the 
DSE calculation should be augmented by the inclusion of pion initial and final 
state interactions. Omitting them, the results described herein are kindred to those 
obtained at tree-level in chiral perturbation theory. Bearing this in mind, in Table 
I4.5l we c ompa re the DSE p redic tions for tt-tt scattering lengths with fits to data using 
tree-leve^3DI anc j two-loorJ ^^ chiral perturbation theory. The comparison shows 
that final state interactions only materially affect the scalar-isoscalar scattering 
length. It is this channel that putatively exhibits a low-mass scalar meson. 

°To complete the proof one uses the axial- vector Ward-Takahashi identity, Eq. 12.81 . the quark-level 
Goldberger-Treiman relation, Eq. 12.421 . and the Bethe-Salpeter equation for the fully-amputated 
quark-antiquark scattering matrix, Eq. (20) of Ref. 1191 . 
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Fig. 4.8. Current-quark-mass-dependence of the isospin-zero and -two scattering amplitudes at 
threshold: j4o(4m^., 0, 0), A2(4m^.,0,0), respectively, calculated by summing the contributions in 
Fig- IHZI (Adapted from Ref. [T3§].) 



It is an opportune point to reemphasise a key aspect of the framework, namely, 
because it works explicitly from a single enunciated kernel, final state interactions 
of the type we have just described can be incorporated cleanly and without over- 
counting. For instance, in Sec. 4.3 we saw that the dressed-quark-photon vertex 
exhibits a real-axis pole when s = —q 2 coincides with a bound state mass and that 
a width was acquired only after including meson rescattering effects. In precisely 
the same manner, the sum of terms in Fig. 14.71 produces the resonance poles one 
expects in 7T-7 T scattering; e.g., a p-meson pole appears in the isospin-one scattering 
amplitude!^ The systematic nature of the truncation means that one can subse- 
quently add loops that represent initial and/or final state interactions between the 
DSE-described rainbow-ladder pions and be certain that these contributions are 
truly new and previously unaccounted for, and just those terms needed to provide 
the strong width of the p-meson. 

It is important to note that in a Ward-Takahashi identity preserving implemen- 
tation of the rainbow-ladder truncation, contributions analogous to the additional 
diagrams in Fig. 14.71 which we have denoted D-r s (s,t,u), Dj- t (s,t,u), can only ap- 
pear in the calculation of amplitudes describing processes with four or more external 



Table 4.5. Rainbow-ladder DSE prediction for it-it scattering lengths, " compared with those 
fitted in analyses of data using chiral perturbation theory at tree-level 1 ^ ! a nd at second-order!!^^ 





°0 


a\ 


a 2 


DSE calc. 


0.17 


0.036 


-0.045 


Tree-level ChPT 


0.15 


0.036 


-0.045 


2nd order 


0.22 


0.038 


-0.044 



DSEs: A tool for Hadron Physics 49 



particles. In those involving only two or three external lines, the addition of a single 
gluon rung connecting any two quarks can be absorbed into the ladder sum that 
generated the vertex or bound state amplitude which joins the two quark lines; i.e., 
it adds to what was already present, and is therefore overcounting. Hitherto, ir-ir 
scattering is the only amplitude with four external lines that has been fully explored 
and many other processes are worthy of study or reexamination. Of particular in- 
terest, perhaps, are the anomalous ^ymm amplitude, wherein the interplay betw een 
the vector and axial- vector Ward-Takahashi identities is known to be i mpor tant!^ 
and, indeed, the Wess-Zumino five-pseudoscalar interaction term itself! I 



4.7. Summary 

We reiterate that every DSE result reviewed in this section is a prediction. All fol- 
low from the strict implementation of the renormalisation-group- improved rainbow- 
ladder truncation defined in Eqs. (|2.12|1 . (|2.13|1 . with the model of the effective cou- 
pling specified in Eq. H^.ljl whose single parameter takes the value in Eq. (13.311 . 
which was chosen to fit /„•, rn 7!l fx, iuk- The analyses demonstrate that the DSEs 
provide a mature approach that links hadron physics experiments with elementary 
properties of QCD and allows for their direct interpretation as critical data on the 
nature of quark confinement and the quark/gluon wave functions of hadrons. The 
success is firm evidence that the quark-quark interaction possesses significant inte- 
grated strength on the infrared domain k 2 < 2GeV 2 , and lattice-QCD simulations 
and DSE studies of QCD's gauge sector are now searching for its origin. 



5. On Baryons 

Contemporary experimental facilities employ large momentum transfer reactions to 
probe the structure of hadrons and thereby attempt to elucidate the role played by 
quarks and gluons in building them. Since the proton is a readily accessible target its 
properties have been studied most extensively and hence an understanding of a large 
fraction of the available data requires a Poincare covariant theoretical description 
of the nucleon. 

The material we have reviewed thus far demonstrates that the properties of 
light mesons are well described by a Poincare covariant rainbow-ladder truncation 
of QCD's DSEs. An extension to baryons begins with a Poincare covariant Faddeev 
equation. That, too, requires an assumption about the interaction between quarks. 
An analysi J ' of the Global Colour ModePm suggests that the nucleon can 
be viewed as a quark-diquark composite. Pursuing that picture yield a Faddeev 
equation, in which two quarks are always correlated as a colour-antitriplet diquark 
quasiparticle (because ladder-like gluon exchange is attractive in the 3 C quark-quark 
scattering channel) and binding in the nucleon is effected by the iterated exchange 
of roles between the dormant and diquark-participant quarks. 

A first numerical study of the Faddeev equation for the nucleon was reported 
in Ref. |144| . and there have subsequently been numerous more extensive analyses; 
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e.g., Refs. |145I146| . In particular, the formulation of Ref. |146j employs confined 
quarks, and confined pointlike-scalar and -axial-vector diquark correlations, to ob- 
tain a spectrum of octet and decuplet baryons in which the rms-deviation between 
the calculated mass and experiment is only 2%. The mod el also reproduces nucleon 
form factors over a large range of momentum transfer J14ZI anc j its de scriptive success 
in that application is typical of such Poincare covariant treatments I W l + 4g | 1 4 ^ | 1 W l 
However, these early successes were achieved without considering the role played 
by light pseudoscalar mes ons. In the context of spectroscopy, studies using the 
Cloudy Bag Model (CBM;P1 

indicate that the dressed-nucleon's mass receives a 
ne gative con tribution of as much as 300-400 MeV from pion self-energy corrections; 
i, e . jl5J, l l52 ] SM + = -300 to -400 MeV. Furthermore, a perturbative study, us- 
ing the Faddeev equation, of the mass shift induced by point like 7 r-exchange be- 
tween quark and diquark constituents of the nucleon obtain^ 5M+ = -150 
to —300 MeV. Such corrections much diminish the value of the 2% spectroscopic 
accuracy obtained using only quark and diquark degrees of freedom. The size 
and qualitative impact of contributions from light-pseudoscalars to baryon masses 
may therefore provide material constraints on the develop ment of a realistic quark- 
diquark picture. This has recently been explored in detaiP^I and we now review 
the findings. 

5.1. Diquarks and a Faddeev equation 

The rainbow-ladder DSE truncation yields asymptotic diquark states in the strong 
interaction spectrum! '*' ™ ' Such states are not observed. Their appearance is an 
artefact of this lowest-order truncation. Higher order terms in the quark-quark 
scattering kerneP (crossed-box and vertex corrections) act to ensure that QCD's 
quark-quark scattering matrix does n ot exhibit singularities that correspond to 
asymptotic diquark bound states P2l^sSI Nevertheless , studies with kernels that don't 
produce diquark bound states, do support a physical interpretation of the masses 
obtained using the rainbow-ladder truncation: m qq plays the role of a confined- 
quasiparticle mass in the sense that l qq = 1 /m qq may be interpreted as a range over 
which the diquark can propagate inside a baryon. These observations motivate the 
following Ansatz for the quark-quark scattering matrix: 

[M qq (k,q-K)l u s = ]T T jP (q;K)A JP (K)T jP (k;-K), (5.1) 
j-p=o+,i+,... 

wherein A J (K) act as diquark propagators and T J (q; K) are Bethe-Salpeter-like 
amplitudes describing the relative momentum correlation of the quarks constituting 
the diquark. 

The validity of Eq. Ij5.1|) is key to the derivation of a quark-diquark Faddeev 
equation for baryons. The simplification is amplified when the summation can be 

P NB. We saw in Sec. 2.2 that such contributions to the quark- antiquark scattering kernel do not 
materially affect the properties of vector and flavour nonsinglet pseudoscalar mesons, a result 
which underlies the successes reviewed in Sec. 4. 
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Table 5.1. Diquark pseudoparticle masses (in GeV) calculated in Ref. 1321 . The magnitudes and 
ordering are ch aracteristic and model independent: see, e.g., Refs. 13411551 and recent lattice-QCD 
estimatesdH (Adapted from Ref. HSH -I 





(ud)+ 


(us)J 




(us)t 


(ss)t 


(uu) 1 


(us) 1 


(ss) 1 


m qq 


0.74 


0.88 


0.95 


1.05 


1.13 


1.47 


1.53 


1.64 



truncated after only a few terms. That is indeed possible, as c an b e argued from 
Table 15.11 which lists calculated diquark pseudoparticle masses. It is apparent 
that for octet and decuplet baryons it should be a good approximation to retain 
only the scalar and axial- vector diquark correlations because they alone have masses 
less than the bound states they will constitute. (Naturally, spin-3/2 baryons cannot 
be described unless pseudovector correlations are retained.) Capitalising on this, 
one arrives at a remarkably simple matrix-integral equation. 

For example, one represents the bound state nucleon by a Faddeev amplitude: 

* = *i + f 2 + ^3 , (5.2) 

in which the subscript identifies the dormant quark and, e.g., ^1,2 are obtained 
from ^3 by a correlated cyclic permutation of all the quark labels. According to 
the above assumption, the individual sub-amplitudes are written as just a sum of 
scalar and pseudovector diquark correlations: 

*3(Pi,a i) r i ) = *§ + +*i + , (5.3) 

with (pi,cti,Ti) the momentum, spin and isospin labels of the quarks constituting 
the nucleon, and P = p\ + p 2 + P3 the nucleon's total momentum. The scalar 
diquark component in Eq. I|5.3|l is 

*§ + (p i ,a i ,r i ) = [r 0+ (\p [12] ;K)}i\z A ° + ( K ) W;PMP)]Z > (5-4) 

where: the spinor satisfies (£7 • P + M)u(P) = = u{P) (17 • P + M), with M 
the mass obtained in solving the Faddeev equation, and is also a spinor in isospin 
space, with ip + = col(l,0) for the proton and ip_ = col(0, 1) for the neutron; and 
K = pi + P2 =■ P{i2}, P[i2] = Pi - P2-, t ■= (—£{12} + 2 P3)/3. The pseudovector 
component is 

* 1+ 7i) = [t'Tf^p^K)]™ All(K) [Al(£;P)u(P)V a l , (5.5) 

where the symmetric isospin-triplet matrices are: t + = (r°+r 3 ) , t° = t 1 , t~ = 
-^(r — r 3 ) , with — <Jy and r 1,3 the usual Pauli matrices. The colour 

antisymmetry of ^3 is implicit in T J . 

The Faddeev equation satisfied by ^3 now reduces to a set of coupled equations 
for the matrix valued functions S, A l u in Eqs. (|5.4I) . 1)5.511 : 

^" r S{£;P)u(P) 



S(k;P)u(P) 
4(fc;P)u(P) 



Al(£;P)u(P) 



(5.6) 
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wherein the kernel is 

M(k,£;P) = 



Moo {M ox )i 
(M 10 % (M n )% 



(5.7) 



with, e.g., 

Moo = T 0+ {k q -£ qq /2; l qq ) S T (£ qq ~k q ) T 0+ (£ q -k qq /2; -k qq ) S(£ q ) A 0+ (£ qq ) , (5.8) 

l q = I + P/3, k q = k + P/3, £ qq = -£+ 2P/3, = -k + 2P/3, and 5 is the prop- 
agator of the dormant quark constituent of the nucleon. The other entries are also 
expressed merely in terms of the Bethe-Salpeter-like amplitudes and propagators. 

It is implicit in Eqs. (|5.6() - (|5.8|l that u(P) is a normalised average of ip± so that, 
e.g., the equation for the proton is obtained by projection on the left with (p\ and 
Mqi generates an isospin coupling between u(P) v+ on the l.h.s. of Eq. (|5.6I) and, 
on the r.h.s., 

V2 At u{P) v - - A , u{P) v+ . (5.9) 

This is merely the Clebsch- Gordon coupling of isospin-1© isospin-^ to total isospin- 
i and means that the scalar diquark amplitude in the proton, (ud) + u, is coupled 
to itself and the linear combination: ^J2 (uu)^ d— (ud)^ u. 

The matrix-valued functions S and A 1 ^ are Bethe-Salpeter-like amplitudes that 
describe the momentum-space correlation between the qu ark a nd diquark in the 
nucleon and can, in general, assume quite complex forms^^ However, reduced 
forms can be used to good effect 

S{£;P) - / 1 (^ > p3)j D + -L( i7 .*_*.pj D ) f 2 (£ 2 ,P 2 ), (5.10) 
4(4 P) = a\{£\P 2 ) lbl ^a\{£\P 2 ) 757^4, (5.11) 

where (Iu)rs = 5 rsi P 2 = —1, i 2 = 1, and, assuming isospin symmetry, aj = a 2 — 
a|, j = 1, 2. In the nucleon's rest frame, /i^ in Eq. I|5.10|l describe the upper, lower 
component of the bound-state nucleon's spinor. It will readily be appreciated that 
a sizable value of /2//1 corresponds to a significant amount of the nucleon's spin 
being stored as quark orbital angular momentum. In this connection the amplitude 
a\ corresponds to a D-w&ve in the nucleon's wave function: in the rest frame, the 
quark and diquark spins are coupled to total-spin |, which must be combined with 
angular momentum L = 2 to obtain a J=| nucleon. The physical content is frame 
invariant because the Faddeev amplitudes are Poincare covariant. 

Equation l|5.6|l can be solved once the kernel is specified and for that Ref. |154j 
adopted the expedient, described in Sec. 3.3, of using algebraic parametrisations of 
the propagators and Bethe-Salpeter amplitudes. The rationale is the same: DSE 
studies of baryons are not yet mature and hence there is merit in employing sound 
simplifications. The dressed-quark propagator of Eqs. (|3.24|I - H3.27|I was used, with- 
out change. The Bethe-Salpeter-like amplitudes describing the relative-momentum 
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Table 5.2. Calculated nucleoli and A massesl ^ ' The results in the first and third rows were 
obtained using scalar and pseudovector diquark correlations: = 0.90 GeV in row 1, m 1 + = 

0.94 GeV in row 3. ("Iq+ = 0.74 GeV, always.) Pseudovector diquarks were omitted in the second 
and fourth rows. u>t are discussed after Eq. 15.161 . and R in and after Eq. 15.191 . All dimensioned 
quantities are in GeV. (Adapted from Ref. 11541 . ) 
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0.45 




1.44 




0.36 


0.35 


2.32 



correlations of quarks within the diquark were represented by 

T 0+ (k;K) = -±-H a Ci l5 iT 2 F(k 2 /u; 2 0+ ), (5.12) 
t' l Tl + (k;K) = ^ri^Cf^H), (5.13) 

with T(y) given after Eq. (|3.26|) . the colour matrix [H C3 ] C1C2 = e clC2C3 and J\f jF , the 
canonical normalisation, fixed by an expression very much like Eq. (|3.34|) : and the 
propagators for the confined diquarks by: 

A° + (K) = -i-^ 2 H 2 + )> (5.14) 
All(K) = (5-15) 

These expressions define a four-parameter model: uijp are diquark masses and ujjp 
are widths characterising the size of the diquark correlation inside baryons. 

The nucleon's Faddeev equation is complete with these definitions. An analo- 
gous equation for the A is readily obtained and does not require additional input. 
The equations can be solved for the nucleon and A masses, and also yield their Fad- 
deev amplitudes with which one can define an impulse approximation to N and A 
elastic and transition form factorsP3SCSZlR e g arc |i n g the four parameters, the scalar 
diquark mass can be taken from Table l5~Tl as may the constraint mj+ /m Q + rs 1.3. 
That leaves the diquark width parameters, ojjp. 

One goal of Rcf. |154j was to illustrate and emulate the success of Ref. |146j 
by showing there are intuitively reasonable values of the parameters for which one 
obtains the nucleon and A masses: Mm = 0.94 GeV, A/a = 1.23 GeV. The results 
of that exercise are listed in Table I5~2*l It is clear that the observed masses are easily 
obtained using solely the dressed-quark and -diquark degrees of freedom described 
above. The first two rows of the table show that the additional quark exchange 
associated with the presence of pseudovector correlations provides considerable at- 
traction. Here it reduces the nucleon's mass by 41% and, of course, the A would 
not be bound in this approach without the 1 + correlation. Furthermore, in agree- 
ment with intuition, the nucleon and A masses increase with increasing tojp . The 
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diquark width parameters are also reasonable. For example, with calculated results 

Z + := 1/cj 0+ = 0.31 fm > l 1+ := l/ui 1+ = 0.17 fin, (5.16) 

these correlations lie within the nucleon (experimentally, the proton's charge ra- 
dius r p = 0.87 fm), a point also emphasised by the scalar diquark's charge radius, 
calculated as described in Ref. |148j : r 2 + = (0.55 fm) 2 . Moreover, defining 
by requiring a least-squares fit of J r (£ 2 /uif 1 2 ) to /i,2(^ 2 ), matched in magnitude at 
i 2 ~ 0, one obtains a scale characterising the quark-diquark separation: 

l q(qq)fi :=l/u h =0.40 fm> 0.15 fm = i Z 0+ . (5.17) 

These scales reveal a significant spatial separation between the dormant quark and 
the diquark-participant quarks while permitting both the quark and diquark to 
remain within the baryon's volume, and thereby provide an intuitively appealing 
picture of confined constituents. For the pseudovector analogue 

Wn =0.36fm> (5.18) 

Finally, the ratio 

R=/ 2 (^ = 0)// 1 (^ 2 = 0) (5.19) 

is a measure of the importance of the lower component of the positive energy nu- 
cleon's spinor. It is not small, a fact that emphasises the need to treat baryons 
using a Poincare covariant framework. As a point of comparison, an analogue in 
the MIT bag model is Rmit := max xe [o,i] ji(2.04ir)/jo(0) = 0.43. 

5.2. Nucleon mass and pion loops 

We have illustrated that an internally consistent and accurate description of the nu- 
cleon and A masses is readily obtained using a Poincare covariant Faddeev equation 
based on confined diquarks and quarks. However, the ttNN and nNA couplings 
are large so it is important to estimate the shift in the masses owing to 7r-dressing. 
That was the primary goal of Ref. |154j but in seeking the estimate a number of 
additional important results were confirmed or established. 

5.2.1. One-loop mass shift 

As a first step, the shift in the mass of a positive-energy nucleon, SM + , was evaluated 
perturbatively; i.e., the effect of just one pion loop on the nucleon's mass was 
calculated. Assuming a pseudoscalar coupling: gNi^^r ■ ttN, g = Mj / ff with M the 
mass of the nucleon in the loop, this can be separated into a sum of three terms: 

5M+ = S F M+ + 5 F M+ + S H M+ , (5.20) 

where 5fM^ represents the Fock diagram with a positive-energy nucleon in the 
loop, 5fM7 is the Z-diagram; i.e., the Fock diagram with a negative-energy nucleon 
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in the loop, and 5hM+ is the tadpole (or Hartree) contribution arising from the 
contact term: [M/{2f%)\ Ntt ■ ttN. 

The loop integrals in Eq. (|5.2()|l must be defined and that was achieved in Ref. 
|154| by implementing a Poincare invariant Pauli-Villars regularisation. It was im- 
mediately apparent that this is equivalent to employing a monopole form factor at 
each irNN vertex: g/(l + k 2 /X 2 ), and since the procedure modifies the pion propa- 
gator it may be interpreted as expressing compositeness of the pion and regularising 
its off-shell contribution. A related effect was identified in Refs. 



Each term in Eq. (|5.20|l is finite after regularisation and straightforward to eval- 
uate numerically. However, some insight can be gained algebraically by supposing 
that the ma ss of t he nucleon in the loop is much greater than any other scale present. 
In that 

6fM+ = -J--i TO ^-m^/+(A)-/ ( +(A) ) (5.21) 

SfM- = -^^ m l(Hml/X 2 ]-l) + mlf w (X) + f- 0) (X), (5.22) 

ShM + = -^K^KAl-lJ + A 3 ), (5-23) 

where all the regularisation-scheme-dependence resides in i)M > 0> which arc 
only functions of the regularisation mass-scale. One sees immediately that, for 
sufficiently small m, > A, 8pM+ + 8pM+ > 0, and it is only when the Hartree 
term is included that 5M + < 0. Furthermore, the leading behaviour of the nucleon's 
mass shift that is nonanalytic in the current-quark mass, m oc \/rn v , is only given by 
Eq. (|5.21l) because DCSB constrains the couplings and ensures a precise cancellation 
between the mm^-terms in Eqs. I|5.22[l . (|5.23() . In general, with the pseudoscalar 
coupling, 6AI + (X) < and decreases monotonically from zero with increasin g A o nly 
because of destructive interference between the tadpole- and Z-diagrams! ' At 
realistic values of the actual value of SM + (A) is determined by the regularisation- 
scale-dependent terms, as visible in Fig. 15.11 

There is no tadpole cont ribut ion to the nucleon's mass shift if one elects to use 
an axial-vector irN coupling! ^ S3 [1/(2/^)] N'y^'y^ f ■ d^ifN. In this case S A M + — 
5 A M+ -\-S A M7 ; i.e., there are two contributions, one with a positive-energy nucleon 
in the loop and the other with a negative-energy nucleon. Some insight is again 
obtained by considering the large-M limit, whereirfi^l 

S A M+ = - ^j^rnl + ml f+ A) (X) + f+ A) (X); (5.24) 

i.e., the same contribution, nonanalytic in the current-quark mass, as in Eq. I|5.21|) . 
but with different regularisation-dependent terms, and 6 M7 oc 1/M because the 
pseudovector coupling suppresses the Z-diagram. The leading nonanalytic contri- 
bution to the nucleon's mass is now unambiguous: it is given by the positive-energy 
nucleon Fock diagram, irrespective of the coupling's character. Implementing the 
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Fig. 5.1. Solid line: Shift in a positive-energy nucleon's mass due to the 0(g 2 ) ^-contribution to 
the self energy. (M = 0.94 GeV, = 0.14 GeV.) SM+{\ = 0.6 GeV) = -0.15 GeV. Dashed line: 
S F M+ Eq. l5"2Tt : dot-dashed line: S A M+ Eq. |5"321 . (Adapted from Ref. 



Pauli-Villars regularisation, a numerical evaluation of 5 M + is straightforward with 
the result depicted in Fig. 15.11 It is evident that 5 A M+ ^ 5fM+, which illustrates 
the difference between the regularisation-dependent terms in Eqs. I|5.21[l and l|5.24|) . 
In addition, although it may not be immediately obvious, 

S A M+ = 5M+ , (5.25) 

which is why there is only one solid curve in the figure. In this outcome one 
has a quantitative verification of the on-shell equivalence of the pseudoscalar and 
pseudovector interactions, in perturbation theory, as long as the interactions are 
treated in a manner consistent with chiral symmetry! 



5.2.2. Mass shift, nonperturbatively 

The effect of infinitely many pion loops can be calculated using a DSE for the 
nucleon's self energy: 

S(P) = 3 J -0j g 2 PV (P, k) A^(P ~ kf) 1 ■ (P - fc) 75 G(k) 7 ■ (P - fc) 75 , (5.26) 

with G-^k) = «7 • k + M + E(P) = i7 • kA(k 2 ) + M + B(k 2 ), where M is the 
nucleon's bare mass, which is obtained, e.g., by solving the Faddeev equation. In 
Eq. H5.2(jfl . A^fc 2 ) = l/[fc 2 + m 2 ] is the pion propagator, and gpy(P,k) is a form 
factor that must describe the composite nature of both the pion and the nucleon. 
The self-consistent solution of Eq. (|5.2ti|) yields A(k 2 ) and B(k 2 ), and therefrom the 
nonperturbative mass shift. 

The ttN vertex function can be calculated using a Poincare covariant model of 
the nucleon, however, a calculation of the mass shift may again be expedited by 
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employing an algebraic parametrisation that is constrained by such studies; e.g.j ^^ 
9pv(P, *) = ^ 9n((P - k) 2 ) g N (P 2 ) g N (k 2 ) , (5.27) 

g 7T (x) = e~ x / A2 , 3at(x) = e~( x+M2 )/ A ™. This model provides for suppression of 
the coupling when either or both the nucleon and pion are off-shell, and in this it 
represents the compositeness of both. The exponential form facilitates an algebraic 
evaluation of many necessary integrals and each term in the product is phenomeno- 
logically equivalent to a monopole form factor 1/[1 + x/X 2 ] if the mass scales are 
related via A = ^2 A. Another advantage of this algebraic form is that it enables an 
elucidation of the precise equivalence between the Minkowski and Euclidean space 
calculations of the mass shift. 

A key aspect of a nonperturbative evaluation of 5M + is that the position of the 
pole in the nucleon's propagator is not known a priori: locating it is the goal, and 
this precludes an algebraic evaluation of the energy-integral that was straightfor- 
ward in the one-loop calculation. In this case one must proceed by first evaluating 
the angular integrals in Eq. Ij5.26|l . which are independent of G(k). That can be 
illustrated with the kernel of the equation for B 



!C B (P 2 ,k 2 ) = Jdn k g 2 PV {P,k) 



2mi 



(P - k) 2 + ml 



(5.28) 



with dQk the usual angular measure. It is evident that /Cg can be considered as 
a sum of two terms. The first is proportional to the angular average of g PV (P, k), 
and using Eq. (|5.27|l that integral can be evaluated exactly: 

f PV (P 2 ,k 2 ) := J dil k g 2 PV ((P~k) 2 ) 

q 2 2(p 2 +k 2 ) n _ o n o A 2 
- e 



n/J - rgliP^glik 2 ) — !^), (5.29) 

where Ii(x) is a modified Bessel function and P — Vp^, k = \/J? . The second term 
is proportional to 

which, in general, cannot be expressed as a finite sum of known functions. However, 
if gpy is regular at P = k and its analytic structure is not a key influence on the 
solution, then the approximation 

uAP\k 2 ) » -^-gi(\P^k 2 \) 9 2 N (p 2 ) 9 2 N (k 2 ) Un k 



2M' WI l;yNK J (P — k) 2 + m 2 



„2 



1 



= ^H9l{\P 2 -k 2 \)9 2 N {P 2 )g 2 N {k 2 ) 2 (5.31) 

--: ~g 2 PV (P 2 ,k 2 ) 7^=1, (5.32) 
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where a — P 2 + k 2 + m 2 , b = 2Pk, is a reliable toolP^ll These preconditions are 
obviously satisfied in this application because the dominant physical effect in nN 
physics is the pion pole and that appears at a mass-scale much lower than those 
present in gpv- Qualitatively identical considerations apply to Kj^. 

The nucleon's mass appears at P 2 < and hence to complete the specification of 
Eq. Ij5.26|) one must define the continuation of the kernels into the timelike region. 
The kernels' primary nonanalyticity is a square-root branch point associated with 
the simple pole in the pion propagator, and in continuing to P 2 < it is necessary to 
include the discontinui ty ac ross the associated cut. That must not be forgotten and 
is readily accomplishe d " ' so that the nucleon's DSE is expressed by two coupled 
integral equations, which we illustrate with that for B 



3 f° 
16tH J Xb 



yA 2 (y)+B 2 (y) 



dyyie B {x,y) 2 , , , (5.33) 



16^ 2 J yy DV yA 2 (y) + B 2 (y) ' 
where Xb = — (\/— x — m^) 2 is the location of the branch point and 

AfC B (x,y) = m n , (5.34) 

xy 

2m 2 

K^B(x,y) = g 2 PV (x,y) - g 2 PV (x,y) ff (5.35) 

a + Vffl — o z 

NB. The A/C^b terms contribute only for P 2 + m 2 < 0. 

With the solutions for A, B in hand, the fully-dressed nucleon mass, Mn, is 
obtained by solving 

M 2 D A 2 (-Ml) = [M + B{-M 2 D )} 2 (5.36) 



and the nonperturbative mass shift is 5M+ = Mb — M . The width s A. , Aw in E q. 
(|5.27(l are constrained by model and phenomenological analyses! f 



A ~ 0.9 GeV , Ajv/A ~ 1.5 - 2.0 (5.37) 
and completing the calculation, 9 Ref. |154j reports a 7rA^-loop induced mass shift 

- 5M+ « (60 - 100) MeV . (5.38) 



Allowing for model-dependence it is therefore safe to say that the it TV-loop reduces 
the nucleon's mass by ~ 10-20%. Extant calculation J ' ' * show that the 
contribution from the analogous 7rA-loop is of the same sign and no greater in mag- 
nitude so that the total reduction is < 20-40%. These same calculations indicate 
that the A mass is also reduced by ir loops but by a smaller amount (~ 50 - 100 MeV 
less). 

9 Three "pions in the air" are sufficient to yield a self-consistent solution and one pion alone provides 
95% of the mass shift. 
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5.2.3. Pions, quarks and diquarks 

The question now is does this contribution materially affect the quark-diquark pic- 
ture of baryons? That may be addressed by solving the Faddeev equations again 
but this time requiring that the quark-diquark component yield higher masses for 
the N and A: M N = 0.94 + 0.2 = 1.14 GeV, M A = 1.232 + 0.1 = 1.332 GeV. 
The resultsPl 

of that exercise are presented in the third and fourth rows of 
Table 15.21 and establish that the effects are material but not disruptive. In this 
case omitting the axial- vector diquark yields Mjv = 1.44 GeV, which signals a 10% 
increase in the importance of the scalar-diquark component of the nucleon. (It is an 
increase because this component now requires less correction. The scalar diquark's 
charge radius was found to be r + = 0.63 fm; i.e., 15% larger.) It also reveals a 
reduction in the role played by axial-vector diquark correlations in the nucleon, since 
now restoring them only reduces the nucleon's core mass by 21%, with it self-energy 
corrections providing the remaining 14%. 

Requiring an exact fit to the N and A masses using only quark and diquark 
degrees of freedom therefore leads to an overestimate of the role played by axial- 
vector diquark correlations: it forces the 1 + diquark to mimic, in part, the effect 
of pions since they both act to reduce the mass cf. that of a quark+scalar-diquark 
baryon. An accurate picture would represent the nucleon as ~ 60% quark+scalar- 
diquark, ~ 20% quark+axial-vector-diquark and ~ 20% pion cloud (irN + ttA). 
This result can significantly impact upon the calculation of quantities such as the 
neutron's charge form factor and the ratio F^q )/Ff (q 2 ), and consequently existing 
Faddeev-equation-based calculations of form factors should be revisited. 



5.3. Nucleon form factors 

The Faddeev equation also pr ovides a bound state amplitude, which is a key element 
in the impulse appr oximat ion^ ^ l ^ ^ used in many studies of the nucleon's electro- 
magnetic and strong form factors! ' U 4 s? l -I4y | l W l However, those studies have all 
overlooked the pion cloud's contribution, which is certainly important at small q 2 
but may also implicitly affect results at larger q 2 by causing some elements in the 
calculations to be overweighted in order to mimic the pion's effect. The studies 
must be updated. Nevertheless, with these remarks in mind, we briefly review a 
topical result. 

References jl48H49j employed a product Ansatz for the nucleon's Faddeev am- 
plitude that retains only the scalar diquark component, *S?® + in Eq. I|5.3|l . This 
exploratory model involved three parameters, which were determined in a least- 
squares fit to G p E (q 2 ) on < q 2 < 3 GeV 2 , and made predictions for a variety of 
other couplings and form factors. Its description of the neutron's charge form factor 
is poor, as it should be given the omission of the axial-vector diquark correlation 
and pion cloud. However, this is the only major quantitative defect; e.g., G p ^(q 2 ) 
are well described. In fact they are probably too well reproduced since 7r-loops 
are expected to contribute ~ 20% to the nucleons' magnetic moments and charge 
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radii.El 

Hence, Refs. |148I149| provide a Poincare covariant DSE model of the nu- 
cleon with a systematic error of < 30%. The model has an important additional 
feature, absent in other approaches: calculated quantities evolve smoothly to their 
perturbative-QCD limit because the DSEs reproduce perturbation theory at weak 
coupling. We saw this explicitly for u*(x) in Sec. 3.3 and F^(q 2 ) in Sec. 4.1. 

We plot results for F%(q 2 )/[KF? (q 2 )], k = F 2 (0), and iiG p E {q 2 )/G p M {q 2 ) in Fig. 
15.21 It was remarked in Ref. I149j that the model agreed semiquantitatively with 
the then-current JL ab d at d^^ but this is the first illustration! ^^*-* ! The data has 
excited some interesl^^ but we judge that here the situation is just as with F 7r (q 2 ) 
[Fig. 14.2] , On the domain hitherto explored the form factors are evolving through 
the region on which infrared phenomena, such as the length-scales defined by bound 
state amplitudes and strongly-dressed quark propagators, are dominant. Thereupon 
a quantitative agreeme nt w ith data is sensitive to model details. Eventually, the 
perturbative behaviour! -*^ ! q 2 F p (q 2 )/ F p (q 2 ) — constant, will become evident but 
that is unlikely until significantly larger q 2 . For instance, with F^i^q 2 ) the perturba- 
tive behaviour is not unambiguously evident until q 2 > 15GeV 2 . It is a challenge 
for quantum field theoretical DSE models of the nucleon, with their unique capacity 
for interpolating between the soft and hard domains, to locate the onset of pertur- 
bative behaviour in nucleon form factors. Given the experience with other elastic 
form factors it would have been a surprise to see that at q 2 < 10 GeV 2 . 
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6. Epilogue 

We have provided a perspective on the contemporary application of Dyson-Schwin- 
ger equations (DSEs) to hadron physics. The keystone of this approach's success 
has always been an appreciation and expression of the momentum-dependence of 
dressed-parton propagators at infrared length-scales. That is responsible for the 
magnitude of constituent-quark and -gluon masses, and the length-scale character- 
ising confinement in bound states. It is now recognised as a fact. Modern hadron 
physics experiments are probing a domain on which this phenomenon underpins ob- 
servable behaviour. The next generation is likely to advance to a region wherein the 
transition to perturbative behaviour takes place. DSE methods, with their unique 
capacity to connect phenomena dominated by soft scales with their perturbative 
limits, will necessarily become increasingly valuable. 

In recent years it has become clear why the simple rainbow-ladder DSE trunca- 
tion has been successful for light vector and flavour nonsinglet pseudoscalar mesons. 
It is the first term in a systematic and nonperturbative scheme that preserves all 
the Ward-Takahashi identities that express conservation laws at an hadronic level. 
Studies in these channels showed that resumming subclasses of diagrams to infinite 
order provides a correction to bound state masses of < 10%. Indeed, with just 
the first order correction to the rainbow-ladder kernel, the calculated masses are 
accurate to 99%. This analysis also explains why the truncation should, and does, 
fail for scalar mesons and points the way toward a systematically improved hadron 
phenomenology 

The existence of a systematic, nonperturbative, symmetry preserving truncation 
scheme has enabled the proof of exact results in QCD. It provides for a straight- 
forward explanation of the dichotomy of the pion as both a Goldstone mode and a 
bound state of effectively very massive quarks. In arriving at that understanding, a 
mass formula for flavour nonsinglet pseudoscalar mesons was exposed which unifies 
the light- and heavy-quark regimes of QCD and provides a qualitative understand- 
ing of lattice simulations and their extrapolation to the chiral limit. 

There have been numerous applications of well-constrained DSE models to 
hadronic phenomena. Among them, a calculation of the pion's valence-quark mo- 
mentum-fraction probability distribution, (x) , highlights the framework's ability 
to provide a description that unifies the soft and hard domains of a given phe- 
nomenon. In this case that has been crucial in exposing a serious discrepancy 
between theory and experiment. In agreement with perturbative QCD, the DSE 
study predicts u*{x) oc (1 — x) 2 in the valence-quark domain. However, extant 
experiments are consistent with u^(x) cx (1 — x). That is profoundly disturbing 
because a verification of the experimental result would even challenge the assumed 
vector-exchange nature of the force underlying the strong interaction. 

The widespread success of a renormalisation-group-improved rainbow-ladder 
model is certainly one of the most significant achievements of the last five years. 
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The model is defined by one parameter, which is an analogue for light-quarks of the 
string tension. That parameter and the two light-quark masses m u = m<2, m s are 
fitted to TOtt, /tt, m,K, fx and everything calculated subsequently is a parameter- 
free prediction. These ab initio calculations have provided a unified description 
of many phenomena, among them: the spectrum of light-vector mesons; tt and K 
electroweak form factors; vector meson transition form factors; and even tt-tt scat- 
tering. And this has not merely been a quantitative success. The studies have 
provided important information on how and where the predictions of perturbative 
QCD become apparent in exclusive processes and, for the pion, tied that to fea- 
tures of low-energy tt-tt scattering and current-algebra constraints. This body of 
work is unique in providing a systematically improvable, nonperturbative, Poincare 
covariant, symmetry preserving approach to hadron physics. 

A significant challenge is to emulate this success with baryons. Studies con- 
ducted hitherto, while preserving these important features, have been exploratory. 
They can and must be improved; e.g., by properly including 7r-cloud effects. To un- 
derstand the data obtained at current and future hadron physics facilities a Poincare 
covariant framework for baryons that naturally expresses the transition from the 
nonperturbative to perturbative domains will be necessary. The early studies sug- 
gest that a Faddeev equation built upon the DSE framework can fill that need. 
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